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Abstract 


The advective transport of atmospheric water vapor and its role in global 
hydrology and the water balance of continental regions are discussed and 
explored. The data set consists of ten years of global wind and humidity 

observations interpolated onto a regular grid by objective analysis. 
Atmospheric water vapor fluxes across the boundaries of selected continental 
regions are displayed graphically. 

The water vapor flux data are used to investigate the sources of 
continental precipitation. 

The total amount of water that precipitates on large continental regions 
is supplied by two mechanisms: (1) advection from surrounding areas external 

to the region and (2) evaporation and transpiration from the land surface 
within the region. The latter supply mechanism is tantamount to the 
recycling of precipitation over the continental area. The degree to which 
regional precipitation is supplied by recycled moisture is a potentially 
significant climate feedback mechanism and land surface— atmosphere interaction, 
which may contribute to the persistence and intensification of droughts. A 
simplified model of the atmospheric moisture over continents and simultaneous 
estimates of regional precipitation are employed to estimate, for several large 
continental regions, the fraction of precipitation that is locally derived. 

In a separate, but related, study estimates of ocean-to-land water vapor 
transport are used to parameterize an existing simple climate model, containing 
both land and ocean surfaces, that is intended to mimic the dynamics of 
continental climates. 
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The horizontal divergence of the vertically-integrated water vapor 
flux, calculated from the original GFDL analysis. Units: 10 cm 

yr" 1 . [Source: Peixoto and Oort (1983)] 

The vertical integral of moisture flux divergence calculated from 
the modified GFDL data. Units: W m* 2 . [Source: Savijarvi 

(1988)] 

Map for each grid point of the lowest level used in vertical 
integration. A value of 1 indicates that the Earth’s surface lies 
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Continental regions selected for graphical analysis of the 
boundary flux of water vapor. The small crosses indicate the 
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South America. Values shown are the average of monthly 
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The total moisture flux (mean motion plus transient eddies) 

across segments of horizontal regional boundaries, for North and 
South America. Values shown are the average of monthly 

average values for June, July and August, 1963—73. 

The total moisture flux (mean motion plus transient eddies) 

across segments of horizontal regional boundaries, for North and 
South America. Values shown are the average of monthly 

average values for September, October and November, 1963-73. 

The moisture flux by transient eddies across segments of 
horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
December, January and February, 1963-73. 
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Figure 2.1: 
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The moisture flux by transient eddies across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 

March, April and May, 1963-73. 

The moisture flux by transient eddies across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for June, 
July and August, 1963-73. 

The moisture flux by transient eddies across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
September, October and November, 1963-73. 

The moisture flux by monthly mean motion across segments of 
horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
December, January and February, 1963-1973. 

The moisture flux by monthly mean motion across segments of 
horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
March, April and May, 1963-1973. 

The moisture flux by monthly mean motion across segments of 
horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for June, 
July and August, 1963-1973. 

The moisture flux by monthly mean motion across segments of 
horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
September, October and November, 1963-1973. 

Simplified model of the atmospheric moisture fluxes over a land 
region. P m and are precipitation of local evaporative and 

advective origin, respectively. 

Fractions of advected and locally-evaporated moisture in the air 
as it moves across a region of length l. 

Plan view of a region, showing segments of the boundary (7) 
across which the atmospheric moisture flux is inward (7.J and 

outward (7 .). 

v 'out 7 


Plan view of a region in which F + , the net moisture influx, is 
not strictly analogous to wu in Budyko’s model. 
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Figure 2.11: 

Figure 2.12: 
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Regions selected for precipitation recycling analysis. The analysis 
technique dictated a rectangular shape; otherwise, the regions 
were chosen on the basis of hydrological and historical interest. 

Eurasian Region. The ratio of locally-evaporated to total 
precipitation, as estimated by the Modified 1-D and the 2-D 
method. 

Eurasian Region. Terms used in estimating P m /P by two 

methods. Atmospheric moisture convergence (-V*<3), Precipitation 
(P) and Evaporation (E) are common to both methods; wu eff /^ eff 

is used in the Modified 1— D method, F+/A in the 2-D method. 

Eurasian Region. Bar graphs of the moisture flux across 

segments of the boundary. Each bar represents the integral of 

the normal component of <3 along the segment at the base of the 
bar. 

North American Region. The ratio of locally-evaporated to total 
precipitation, as estimated by the Modified 1-D and the 2-D 
method. 

North American Region. Terms used in estimating P m /P by two 

methods. Atmospheric moisture convergence (-V*<3), Precipitation 
(P) and Evaporation (E) are common to both methods; wu eff /^ ff 

is used in the Modified 1-D method, F + /A in the 2-D method. 

North American Region. Bar graphs of the moisture flux across 
segments of the boundary. Each bar represents the integral of 

the normal component of 3 along the segment at the base of the 
bar. 

South American Region. The ratio of locally-evaporated to total 
precipitation, as estimated by the Modified 1-D and the 2-D 
method. 

South American Region. Terms used in estimating P m /P by two 

methods. Atmospheric moisture convergence (-V*<3)> Precipitation 
(P) and Evaporation (E) are common to both methods; 

is used in the Modified 1-D method, F^/A in the 2— D method. 
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Figure 3.1 
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South American Region. Bar graphs of the moisture flux across 
segments of the boundary. Each bar represents the integral of 

the normal component of <3 along the segment at the base of the 
bar. 

African Region. The ratio of locally-evaporated to total 
precipitation, as estimated by the Modified 1— D and the 2— D 
method. 

African Region. Terms used in estimating P m /P by two 

methods. Atmospheric moisture convergence (-V-<3), Precipitation 
(P) and Evaporation (E) are common to both methods; wu eff /^ eff 

is used in the Modified 1— D method, F"^ / A in the 2— D method. 

African Region. Bar graphs of the moisture flux across segments 
of the boundary. Each bar represents the integral of the normal 

component of (3 along the segment at the base of the bar. 

a) Plan view of a region with convergent moisture flux, and with 
constant f<3|| along the inflow boundaries. By symmetry, the 
average direction of moisture flux (u) is due eastward. 

b) The difference between the total boundary integral of influx 
(F + ) and the integral of boundary influx resolved parallel to u 
t(wu) cff Ay)]. 

Streamlines of the total aerial runoff, <3; each barb represents 2 g 
kg' 1 rn S' 1 : b, DJF; g, JJA; C, year. [Source: Peixoto and Oort, 1983] 

The continent and ocean regions used to estimate the transport 
parameter, x, from observations: Experiment 1. 

The continent and ocean regions used to estimate the transport 
parameter, x, from observations: Experiment 2. 

X and -V-<3 for six different runs of the screening model. For 
each run, all other model parameters were held fixed while x was 

varied; the model was integrated over 250 days for each (-V-<3, 
X) pair. 

Estimates of x and -V-<3 from the observational data for each of 
twelve months. The observations that do not satisfy model 
assumptions are labeled by month. 
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Figure A. 2 
Figure A. 3 


Estimates of x an d -V*($ from the observational data, for selected 
months. The observations that did not satisfy model assumptions 
have been excluded. 

Observational estimates of (-V-3, x) compared with model 

output (-V-<3> x) produced by varying x with other 
parameters fixed, Run a. 

Observational estimates of (-V-4 x) compared with model 

output (-V-3, x) produced by varying x with ^ other 

parameters fixed, Run b. 

Observational estimates of (-V-3, x) compared with model 

output (— V • cl, x) produced by varying x with ^ other 

parameters fixed, Run c. 

Observational estimates of (-V-3, x) compared with model 

output (-V-3, X) produced by varying x with 311 other 

parameters fixed, Run d. 

Observational estimates of (-V-4 x) compared with model 

output (-V-3, x) produced by varying x with all other 

parameters fixed, Run e. 

Observational estimates of (-V-3, x) compared with model 

output (-V-3, x) produced by varying x with all other 

parameters fixed, Run f. 

Region Testl. Bar graphs of the moisture flux across segments 
of the boundary. Each bar represents the integral of the normal 

component of along the segment at the base of the bar. 

Region Test2. Bar graphs of the moisture flux across segments 
of the boundary. Each bar represents the integral of the normal 

component of 3 along the segment at the base of the bar. 

Definition sketch for Section A.l. 

Definition sketch for Section A. 2. 

Definition sketch for Section A. 3. 
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Chapter 1 

Atmospheric Water Vapor Transport 


1.1 Introduction 

Water vapor is a variable constituent of the atmosphere; it accounts for 
0 to 4 percent by volume of tropospheric air. The saturation vapor pressure 
decreases with temperature; therefore the vapor content of the atmosphere 
decreases with height, with the result that atmospheric water vapor is 
concentrated in the lower troposphere and is negligible above the tropopause. 
The atmosphere contains about 12 * 10^ km^ equivalent of liquid water, 
enough to cover the Earth’s surface to a depth of about 2.4 cm (assuming a 
smooth, spherical Earth). The average residence time of a water molecule in 
the atmospheric reservoir is on the order of ten days. 

Water vapor is transported in the atmosphere by molecular and turbulent 
diffusion, convection, and advection. The vertical flux of evaporation from the 
oceans and the land surface into the atmosphere is accomplished by diffusion, 
and the vapor is vertically mixed by diffusion and convection. On the 
temporal and spatial scales of lateral transport, advection is the dominant 
mechanism. On these scales, water vapor can be treated as a scalar 
admixture advected by the horizontal wind. 

The atmospheric water vapor transport vector has been an object of 
study, on both the global and the regional scale, since the availability 
(following World War II) of the aerological data required for vapor flux 
computations. Benton and Estoque (1954) used the vapor flux vector 
divergence field in computations of evapotranspiration over the North American 
continent for the calendar year 1949. In his 1956 Presidential Address to the 
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Royal Meteorological Society, Sutcliffe called for increased scientific inquiry 
into the problem of water balance and the general circulation of the 
atmosphere, and recommended wider application of the techniques introduced 
by Benton and Estoque. Starr and Peixoto (1958) found centers of divergence 
of atmospheric water vapor flux over many of the world’s deserts, and inferred 
from mass balance considerations the existence of net subsurface inflow to 
these regions. Lufkin (1959) made the first estimate of water vapor transport 
from the oceans to the continents as a function of latitude. Hastenrath (1966) 
analyzed water vapor flux over the Central American seas as part of a more 
comprehensive study concerned with the general circulation and energetics in 
that area. Rasmusson (1967, 1968) investigated water vapor flux over North 
America and the Central American Sea as part of a continental water-balance 
study; this work was followed by a regional study of the hydrology of eastern 
North America (Rasmusson, 1971). Rasmusson (1977) recommended the 
application of vapor flux data in the routine computation of regional water 
balances over areas large enough to control the relative error of the resulting 
flux divergence estimates (he recommends 10* km 2 or more). He also provided 
guidelines for operational vapor flux computations within the context of the 
World Weather Watch. Peixoto and Oort (1983) used the divergence field to 
study connections between the atmospheric and surface branches of the 
hydrologic cycle. Chen (1985) analyzed water vapor transport and 
maintenance during the solstice seasons to explore how the high water vapor 
content of certain tropical areas is maintained by the large scale atmospheric 
circulation. 

This thesis explores the role of the advective transport of atmospheric 
water vapor in two different aspects of global hydrology and the water balance 
of continental regions. This Chapter introduces some definitions and concepts 
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that are common to both investigations, and describes the observational data 
used in this study. Chapter Two addresses the question of the sources of 

precipitation water in selected land regions. Some of the water vapor that 
condenses into precipitation over a region is derived from vapor advected from 
the surrounding areas external to the region itself. The remainder is derived 
from surface evaporation within the region. The latter source is governed by 
regional surface hydrologic processes and it constitutes an important mechanism 
for land surface-atmosphere interaction. In Chapter Three, the aerological 
data are used to assess the validity of a parameterization for moisture 
advection in a spatially-averaged climate model. The simplified climate model 
represents the atmospheric branch of the hydrologic cycle over a land surface 
and an ocean region as exchange between two lumped reservoirs. The analysis 
in Chapter Three compares the parameter characterizing this exchange in the 
model with those derived from the aerological data. 

1.2 Basic Concepts 

The Vertically -Integrated Vapor Flux Vector 

The atmosphere’s specific humidity and the zonal and meridional 
components of the wind are variable in both space and time. In the 
geophysical coordinate system, the four-dimensional domain of these variables 
is defined by <p (latitude), A (longitude), p (pressure), and t (time). Following 
Peixoto (1973) and Peixoto and Oort (1983), the dimensionality of the domain 
is reduced to two by defining a vertically-integrated, time-averaged vapor flux 
vector. 
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At a given point in space and time, a vector field of the advective 
transport of water vapor by the horizontal wind is defined: 

qtf(A,^,p,t) = (qu)i + (qv)j (1.1) 

in which q is the specific humidity [g H 2 0 vapor per kg air]; ^ is the 
horizontal wind vector [m s *]; u is the zonal wind component (positive 
eastward); v is the meridional wind component (positive northward); and i 
and j are the unit vectors in the zonal and meridional directions, respectively. 
The horizontal transport vector q\^ has units g kg -1 m s -1 . As a product of 
terms, q\f increases with either velocity or specific humidity; for example, 
slow-moving moist air may transport water vapor at the same rate as less 
moist but fast-moving air. 

For a vertical atmospheric column with a unit base, the 

vertically-integrated horizontal vapor flux, 3, is obtained by taking the 
mass-weighted vertical integral of (1.1): 

p 0 

3(A,&t) = J qV = Q A i + Q^j (1.2) 

0 

in which p Q is the surface pressure; g is the acceleration of gravity; and 
and are, respectively, the zonal (positive eastward) and meridional (positive 
northward) components of C$. The vapor flux vector 3 has units g s -1 m -1 . 
As a vertically-integrated quantity, 3 expresses the magnitude and direction of 
the net transport of water vapor through the depth of the atmosphere above a 
point on the earth’s surface; a small value of either component does not 
necessarily imply negligible transport at all levels in the atmosphere. 

Equation (1.2) may be averaged over a time period r to compute the 
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corresponding mean values and Q^, where the overbar denotes the 

time-average operator: 


Q^(-M) 



dt 


0 


(1.3a) 


and 


Q = \ j dt 

0 


rT 


(1.3b) 


The time-averaged, vertically-integrated moisture flux vector <3 has units 
g s -1 m -1 ; it is a function of two horizontal space variables, latitude and 
longitude. 


The Role of <3 in the Basic Equations of Hydrology 

The mass conservation equation for water vapor in the atmosphere may 
be expressed as 

|y + V-3 = E — P (1.4) 

in which W is the water vapor storage within the atmospheric column; E is 
the rate of evapotranspiration into the base of the column; P is the rate of 
precipitation from the column, also measured at the surface; and is the 
divergence of the lateral vapor flux, with ^ as defined above. W, also known 
as the precipitable water in the column, may be expressed as either water 
vapor mass per unit surface area [kg m -2 ] or equivalent depth of liquid water 
[m]. The use of the latter unit is more common in hydrologic practice; 
correspondingly, P, E, and have units m s' 1 . 

For an atmospheric control volume, bounded by a conceptual closed 
vertical wall and overlying a region of the Earth’s surface having area A, (1.4) 
takes the form 
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(1.5) 


<3T> + <*-<3> = < E - ■>> 

where the angled brackets denote a space average: 


in which h(x) is any function of location. Applying Green’s theorem allows 
(1.5) to be written 

<fr> = < e ■ p > ■ H ^ ■ “7 d7 (L7) 

in which n^ is the outward unit vector normal to the horizontal boundary (7) 
of the region. Equation (1.7) states that the time rate of change in water 
vapor storage within the atmospheric volume is equal to evapotranspiration 
minus precipitation less the net lateral outflow through the vertical boundaries 
of the control volume. 

Equations (1.5) and (1.7) express the mass balance of water in the 

atmospheric branch of the hydrological cycle. An excess of evapotranspiration 
over precipitation is balanced by the change in water content of the 
atmosphere and the net aerial runoff across the boundaries of the region. A 

region in which the mean divergence of < 3 , (V*C$), is positive is a source of 

water vapor to the rest of the atmosphere, while a region with negative 
<V-3>, a region of convergence, is a sink of atmospheric water vapor. 

The mean net transport of water into a region by horizontal atmospheric 
motion is equal to -(V-3) and the total transport into the region is the closed 
path integral - 0 <3-n^ d7, with the normal to the boundary defined outward 
as in (1.7). In analyzing and modeling the atmospheric water balance, the 
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expression "divQ" (omitting the negative sign) is frequently used in reference 
to either of these quantities. 

In the land surface branch of the hydrologic cycle, the mass balance can 
be written in the simplified form, 

<f> = < p - E > + < R i„ - R „„.> ( L8 > 

in which S is the mass of water stored in the surface and subsurface 

reservoirs of the region; R.^ is the rate of water influx across the horizontal 

boundaries of the region by surface and subsurface lateral transport; and R qu( . 
is the rate of water efflux by these lateral processes. If the region’s 
boundaries are surface drainage divides, and if the subsurface drainage 
corresponds to that on the surface, then (R- n ) = 0 an d the term 

(R. - R ) is simply the rate of runoff per unit area; such a simplification 

' in out' 

is also the case when the study region is an entire continent. For smaller 

land regions that are not delineated along natural drainage divides, surface and 
subsurface lateral transport into a region may be significant. 

Equation (1.8) may be rearranged as follows: 

<P - E) = - <E - P> = <R oul - R.J + <g> (19) 

which states that an excess of precipitation over evapotranspiration is balanced 
by net surface and subsurface runoff and/or increased terrestrial storage in the 
region. 

The term (E - P) is common to (1.7) and (1.9) and thus provides the 
link between the terrestrial and atmospheric branches of the hydrologic cycle: 

< R out " R in> + <H> = “ - <!t> • (L10) 
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Therefore, if the aerological data are known, as well as runoff and 
precipitation, (1.10) and (1.7) may be used to estimate the average rate of 
change of groundwater storage and the average evaporation rate, two quantities 
that are difficult to measure. The accuracy of such an estimate is highly 
sensitive to the quality of the measurements. Typically, ^ and ^ are 
smaller than the other terms by at least one order of magnitude, and the 
smaller quantity is estimated from a difference of larger numbers. Accuracy is 
improved by enlarging the path of the boundary integral (1.7) which leads to 
the recommendation of Rasmusson (1977) that the region should have an area 
of at least 10 6 km 2 


Modes of Transport 

Any atmospheric variable, f, may be expressed in terms of a time mean 
and a perturbation from that mean: 

f(A,0,p,t) = f(A,^,p) + f'(A,0,p,t) (1.11) 

in which the overbar denotes the mean value over a specified time interval. 
By definition, the mean value of the perturbations is zero (f = 0). 

For example, expressing the zonal component of the horizontal wind and 
the specific humidity in this manner, the time mean of the zonal component of 

the horizontal vapor transport vector, qV, contains cross-product terms: 


qu(A,0,p) = (q + q')(u + u') 


( 1 . 12 ) 


= q u + q u' + q'u + q'u' . 
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(1.13) 


Because q and u are constants, and u' and q' vanish by definition, 


qu = q u + q'u' . 


The first term on the right-hand side of (1.13) represents the advection of 
mean specific humidity by mean zonal atmospheric motion. The second term 
is a correlation term; for example, q'u' is positive for a given location if 
winds that are more westerly (or less easterly) than the mean wind (positive 
u') are moister than average (positive q'), negative if winds that are more 
easterly (or less westerly) than the mean wind (negative u') tend to be 
moister than average, and zero if no correlation exists between the zonal wind 
and the specific humidity. Likewise, for the meridional component of vapor 
transport, 

qv(A,0,p) = q v + q'v' . (I- 14 ) 

By convention, the terms on the right-hand sides of (1.13) and (1.14) 
are called, respectively, transport by mean motion and by transient eddies. 
The distinction between mean and eddy transport depends upon the averaging 
period. 

The time-averaged vertically-integrated water vapor flux vector, ^ , 
contains a mean and a transient eddy term in each of its components, i.e., 


and 




‘‘"i 2 


+ 




(1.15a) 



0 



o 



(1.15b) 
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Globally, the zonal component of 3 tends to exceed the meridional 
component by an order of magnitude. The mean motion dominates in Q^, 
whereas in Q^, the mean and transient eddy terms are of comparable magni- 
tude. Transient eddy transport is particularly important at mid-latitudes and 
during the active winter season, where the eddies play a critical role in the 
poleward transport of sensible, as well as latent, heat (Peixoto and Oort, 
1983). 

A similar decomposition can be performed in the spatial dimension by 
averaging with respect to longitude around a belt of constant latitude. The 
resulting zonal mean and standing eddy terms are of great interest and value 
in studies of vapor transport over the globe as a whole. The present study 
focuses on the water balances of continental regions on a spatial scale that is 
too small to justify such spatial decomposition. However, the distinction 
between the time mean and transient eddy modes of transport will figure in 

the computation of <3 and the analysis of vapor flux into and out of regions. 

1.3 Observational Data and Computational Method 

Aerological Data 

The humidity and wind data used in this study are from a data set 
provided by the Geophysical Fluid Dynamics Laboratory (GFDL) of NOAA at 
Princeton University through the courtesy of Dr. Abraham Oort. The GFDL 
data set consists of ten years of observations, transformed into a gridded form. 
The measurements were taken during the period May 1963 through April 1973. 
For the years 1963-1968, the observations were once daily, at 00 GMT, 
although a small number of 12 GMT observations were included. For the 
years 1968-1973 both 00 and 12 GMT data were used (Oort, 1983, pp. 5-8). 
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The aerological data were interpolated onto the regular grid by means of an 
objective analysis scheme, using the zonal average of data in a latitudinal belt 
as a first approximation. The objective analysis scheme, known as CRAM 
(Conditional Relaxation Analysis Method) is described by Harris et al. (1966) 
and Rosen et al. (1979). 

Values of q, u, v, q'u' and q'v' are given at each node of a 2.5° 
latitude by 5° longitude grid, and at each of 11 pressure levels; the overbars 
denote monthly averages, with q'u' and q'v' being the horizontal transient 
eddy fluxes in the zonal and meridional directions, respectively. 

A known problem with the GFDL data set is the existence of suspicious 
divergences in the mean winds at mid-latitudes. Savijarvi (1988) used the 
vorticity equation and smoothing techniques to force a mass balance of the 
mean winds in the GFDL data. The effects of Savijarvi’s correction on the 
water vapor flux divergence field are illustrated in Figures 1.1 and 1.2. Figure 
1.1, taken from Peixoto and Oort (1983), shows V.($ for the solstice seasons as 
calculated from the original data. Figure 1.2, from Savijarvi (1988) shows the 
corrected V-<3 for the same seasons. (The figures are not given in the same 
units; multiplying values in Figure 1.1 [10 cm yr _1 ] by 6.3 converts them to 
[W nr 2 ], compatible with Figure 1.2.) Over the continents, several changes 
are noteworthy, including the appearance in Fig. 1.2 of a tongue of 
convergence over Mexico and the U.S. Southwest in JJA, the narrower belt of 
convergence over Africa in JJA, and the disappearance of an irregularly-shaped 
zone of convergence over South America below 20° S. The "unexpected intense 
center of divergence" just west of India in JJA noted by Peixoto and Oort 
(1983) persists in the corrected divergence field. Unfortunately, Savijarvi’s 
improvements to the GFDL data set were brought to the author’s attention 
too late to be incorporated into the present study. 
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Figure 1.1: The horizontal divergence of the vertically-integrated water vapor 

flux, calculated from the original GFDL analysis. Units: 10 cm 

yr _1 . [Source: Peixoto and Oort (1983)] 
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Figure 1.2 The vertical integral of moisture flux divergence calculated from 
the modified GFDL data. Units: W nr 2 . (Source: Savijarvi 

( 1988 )] 
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Alestalo (1983) applied a constant error method to the station radiosonde 
data used in a study of the European water vapor budget. The following is a 
summary of Alestalo’s correction techniques. 

At a given pressure level, the divergence of the time mean water vapor 
transport vector can be expressed as 


V- qtf = q(V-tf) + \^-Vq + V-q'\f' . (1.16) 

The divergence of the time mean wind can be expressed as 


V-tf = V-tf + (V-tf)" 


(1.17) 


in which /\ denotes a vertical mean and ( )" a perturbation from that mean. 
Incorporating (1.17) into (1.16) and integrating vertically gives the result 







q(V-\^)" 


d E 
g 



(1.18) 


Invoking the boundary condition that time mean vertical velocity vanishes at 
the upper and lower boundaries, and assuming no accumulation of air mass in 

the column, must vanish as required by mass conservation. However, this 
is not the case with real observations and analyzed data, due to sampling 
error, variability in the physical quantity being measured, and errors 
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introduced by interpolation and numerical integration. The non-zero 
divergence of the time-mean wind is eliminated by computing a corrected 
moisture flux divergence at each level, as follows: 


V-qtf 


corrected 


= V-qtf 


computed 


- q 


(1.19) 


The correction in (1.19) ensures that air mass is conserved for the column as 
a whole. However, at individual levels the divergence estimates may still be 
erroneous. Since water vapor advection is largely confined to the lower few 
atmospheric layers, the errors at individual levels may still be a problem, 
although the whole column is non-divergent in mass. 

In the present study, no mass-balance corrections were applied to the 
aerological data. The nature of the analysis techniques applied herein would 
require corrections to the vapor flux components, and not only to the vapor 
flux divergence. The development of an appropriate correction method is of 
vital importance, as the results will show. 


Computational Method 

The components of the vertically-averaged water vapor flux vector, as 
defined in (1.15), were evaluated by trapezoidal rule integration. For each 

node, the surface pressure, p , was set equal to its mean annual value, with 
the beginning pressure level for integration selected according to Figure 1.3. 
Figure 1.3, taken from Oort (1983), indicates the annual mean surface pressure 
at the grid points in terms of the pressure levels of the analyzed data. 

A value of 1 indicates that the earth’s surface lies between 1000 and 950 mb, 

a value of 2, between 950 and 900 mb, and so on. Vertical integration was 
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Figure 1.3: Map for each grid point of the lowest level used in vertical 

integration. A value of 1 indicates that the Earth’s surface lies 
between 1000 and 950 mb, a value of 2 between 950 and 900 mb, 
a value of 3 between 900 and 850 mb, a value of 4 between 850 
and 700 mb, and a value of 5 between 700 and 500 mb. 
[Source: Oort (1983)] 
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begun at the next standard pressure level above the surface pressure, i.e., at 
1000 mb for a surface pressure greater than or equal to 1000, at 950 mb for 
surface pressure lying between 1000 and 950 mb, and so on. 

Precipitation Data 

The precipitation data used in this study are from the World Monthly 
Surface Station Climatology (and Associated Datasets) distributed by the 
National Center for Atmospheric Research (NCAR) in Boulder, Colorado. 
Monthly station precipitation data are given for over 3900 different World 
Meteorological Organization stations through 1987, with data for some stations 
going as far back as the mid-1700’s. From this data set, the precipitation 
records for May 1963 through April 1973 were selected, in order to obtain 
precipitation data that are temporally compatible with the GFDL aerological 
data. 

1.4 Boundary Flux Estimation 

The closed-path line integral in (1.7) may be expressed as a sum of 
sub-integrals over segments of the boundary: 

m 

o <3.n 7 d7 = ^ 
k = l 

in which x k and x fc+1 are longitude-latitude points in sequence along the 
boundary, and m is the number of segments. For a closed path, x m+1 must 
be the same point as x . Each sub-integral has units kg s _1 and represents 
the time rate of moisture mass flux across that segment of the boundary. 
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For clarity in discussing the fluxes of moisture over the boundaries of 
regions, additional notation is introduced. The vertically-integrated moisture 
flux vector is decomposed into its mean motion and transient eddy terms, 


^ + ^TE 

where 

3 te (a,*) s f s f- 

0 

and 

S | S q ^ g^ ' 

o 

The boundary sub-integrals in (1.20) are abbreviated as 



( 1 . 21 ) 

(1.22a) 

(1.22b) 

(1.23) 


and, analogous to (122), 


and 


By definition, 




f k = 


1- + f 


M,k 


TE,k 


(1.24a) 


(1.24b) 


(1.25) 


Consistent with the definition of n^ as the outward unit vector, f^ is positive 
for moisture flux out of the region, and negative for flux into the region. 
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If a rectilinear boundary is defined by connecting nodes in the GFDL 
data grid, then each sub-integral, f^, represents the moisture flux across either 
an east-west or a north-south line segment on the Earth’s surface. Such 

rectilinear boundaries are well suited to graphical display of the moisture flux 
data; this technique is now introduced. 

Figure 1.4 shows the rectilinear boundaries of North and South America 
superimposed on a cylindrical equidistant projection of the continents. The 
small crosses indicate the nodes of the GFDL grid (every 2.5° latitude and 5° 
longitude), and the solid lines the boundaries of the study regions. The nodes 
define the endpoints of the boundary segments, over which the sub-integrals in 
(1.24) were evaluated by trapezoidal rule integration for each of the twelve 
months of the year. The results are shown in a set of figures below. In 
these figures, the magnitude of the moisture flux, f^, across segment k of the 
boundary is shown by a bar perpendicular to segment k, with its base at the 
segment. The length of the bar is proportional to f^ by the scale indicated, 
and the bar points into or out of the region according to the direction of net 
moisture flux across segment k. In addition, those bars representing flux into 
the region (negative f^) are cross-hatched to further distinguish them from 
bars representing efflux (f^ positive). The total influx of atmospheric moisture 
into the region is proportional to the total length of cross-hatched bars, and 
the efflux is proportional to the total length of white bars. 

It is important to note that the bars do not represent vector 
components; each represents the integral of the normal component of the flux 
vector over the length of the segment. The graphical scale, in units of kg s* 1 , 
refers to the bars and has no relevance to the continental outlines, which are 
not drawn to scale and are in fact highly distorted by the cylindrical 
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Figure 1.4: Continental regions selected for graphical analysis of the 

boundary flux of water vapor. The small crosses indicate the 
nodes of the GFDL data grid (5° longitude and 2.5° latitude). 
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equidistant projection. Figure 1.5 shows a simple example, illustrating the 
meaning of the bars in the set of figures to follow. Consider a constant Q, 
directed at angle a = tan'^O.S) (so that is twice Q.) and two boundary 
segments, one of 5 degrees longitude, the other 2.5 degrees latitude. At the 
Equator (case a), the east-west boundary segment is twice as long as the 
north-south segment; however, the normal components of the flux vector across 
the boundaries have the inverse relationship; therefore, in this case, f^ is equal 
to f 2 . At 60° North (case b), the two segments have the same length, with 
the result that f is equal to one-half of f 2 - The right-hand side of Figure 
1.5 shows how these results would be plotted in longitude-latitude space 
(equivalent to the cylindrical equidistant projection), in a figure such as 1.6 
through 1.17. 

Figures 1.6 through 1.9 show the total moisture flux (f) across the seg- 
ments of the rectilinear continental regions defined in Fig. 1.4, for the months 
December, January and February (DJF); March, April and May (MAM); June, 
July and August (JJA); and September, October and November (SON). Fig- 
ures 1.10 through 1.13 present the transient eddy moisture flux (f TE ) and 
Figures 1.14 through 1.17, the mean motion flux (f M ) for the same seasonal 
sequence. For Figures 1.14 through 1.17, values of f TE ^ and f M ^ were com- 
puted, then the values for three months were simply averaged. Therefore, the 
transient eddy values plotted in Figures 1.10-1.13 represent moisture transport 
by motions on time scales less than one month during the three month period, 
and the mean motion values plotted in Figures 1.14-1.17 represent transport 
on the time scale of one month. 

The total flux of water vapor onto these continents reflects the 
circulation of the atmosphere and the distribution of temperature at the 
Earth’s surface and in the lower troposphere. The pattern of the equatorial 
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Figure 1.5: A simple example illustrating the meaning of the bars in Figures 

1.6 through 1.17. For constant ($, the magnitude of f depends 

upon the normal component of ($ and the length of the segment. 
The right-hand side shows how these results would appear in 
Figures such as 1.6 through 1.17. 
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Figure 1.6: The total moisture flux (mean motion plus transient eddies) 

across segments of horizontal regional boundaries, for North and 
South America. Values shown are the average of monthly 
average values for December, January and February, 1963-73. 
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Figure 1.7: The total moisture flux (mean motion plus transient eddies) 

across segments of horizontal regional boundaries, for North and 
South America. Values shown are the average of monthly 
average values for March, April and May, 1963-73. 
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1.8: The total moisture flux (mean motion plus transient eddies) 

across segments of horizontal regional boundaries, for North and 
South America. Values shown are the average of monthly 
average values for June, July and August, 1963-73. 
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1 . 9 : The total moisture flux (mean motion plus transient eddies) 

across segments of horizontal regional boundaries, for North and 
South America. Values shown are the average of monthly 
average values for September, October and November, 1963-73. 
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Figure 1.10: The moisture flux by transient eddies across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
December, January and February, 1963-73. 
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Figure 1.11: The moisture flux by transient eddies across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
March, April and May, 1963-73. 
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Figure 1.12: The moisture flux by transient eddies across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for June, 
July and August, 1963-73. 
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Figure 1.13: The moisture flux by transient eddies across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
September, October and November, 1963-73. 
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Figure 1.14: The moisture flux by monthly mean motion across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
December, January and February, 1963-1973. 
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Figure 1.15: The moisture flux by monthly mean motion across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
March, April and May, 1963-1973. 
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Figure 1.16: The moisture flux by monthly mean motion across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for June, 
July and August, 1963-1973. 
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Figure 1.17: The moisture flux by monthly mean motion across segments of 

horizontal regional boundaries, for North and South America. 
Values shown are the average of monthly average values for 
September, October and November, 1963-1973. 
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easterly trade winds and the mid-latitude westerlies is clearly reflected in the 
zonal transport. The warm air moving off the equatorial Atlantic Ocean 
delivers moisture at an extremely high rate over the northeast coast of South 
America. The large magnitudes of the segment fluxes both onto and off 
Patagonia (the southern tip of South America) are attributable to the high 
wind velocities in the "roaring forties" latitudes; at these latitudes the winds 
circle the globe almost entirely over ocean, which provides much less resistance 
to air flow than does land, and at the same time allows development of high 
relative humidity. These bars are longer by half in the Southern Hemisphere 
summer (DJF) than in the winter (JJA), due to the higher saturation vapor 
pressure of the warmer air. On the other hand, the very low landward values 
of f along the west coast of South America between 35 S and 20 S are due to 
the barrier posed by the Andes range. Computationally, as can be seen from 
Figure 1.3, the vertical integration of <3 for these grid nodes begins fairly high 
in the atmosphere, where specific humidity is low. Physically, the small bars 
reflect the loss of water from the air due to adiabatic cooling and resulting 
condensation as the air is forced up over the mountain range. Two major 
streams of atmospheric water vapor onto North America are evident: flow 

from the Pacific Ocean and from the Gulf of Mexico. The latter source is 
strongest in the Northern Hemisphere spring, MAM. 

The meridional transient eddy flux of moisture onto these continents is 
poleward in all seasons (Figures 1.10-1.13). In North America, f TE is most 
prominent off the Gulf of Mexico, where it is greatest in the Northern 
Hemisphere winter, accounting for more than half the total flux across the 
boundary between 100 W and 80 W. Because the rectilinear regional 
boundary does not exactly follow the eastern coastline of North America, it is 
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not clear whether the meridional bars on those boundary segments actually 
correspond to onshore transport of moisture. The apparent unimportance of 
the transient eddy mode for the South American continent is possibly 
attributable to data sparsity. 

Figures 1.14 through 1.17 do not include any new information because f^ 
can be obtained by subtracting f TE from f. However, they are included for 
comparison purposes. In particular, f M from the Gulf of Mexico is quite large 
in MAM and JJA; recalling that f TE is quite small across that boundary in 
JJA, it is interesting to note that the year-round supply of moisture from the 
Gulf of Mexico is accomplished by different modes in different seasons. The 
zonal mean vapor transport off the Gulf of Mexico reverses direction in the 
course of the year, from westerly in DJF to easterly in JJA, an important 
factor in the annual march of precipitation in the Southwest USA (Hastenrath, 
1966). Across a number of boundary segments along the east coast of North 
America, f M and f TE are of opposite sign; this opposition holds mostly for 
meridional transport, although in DJF (Figs. 1.10 and 1.14), a small zonal f TE 
is opposite f M between 4 S and 60 N. In South America, meridional values of 
f TE and f M tend to be of opposite sign along the southeast coast. 

Bar graphs of the boundary flux of ($ are a useful technique for the 
presentation of moisture flux data. Although it would entail substantial 
interpolation between nodes, the ability to define more realistic continental 
boundaries would improve the method. 
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Chapter 2 

Estimation of Continental Precipitation Recycling 


2.1 Introduction 

The total amount of water that precipitates on a region of land comes 
from two sources: 1) water vapor transported into the region by atmospheric 

advection (as described in Chapter One), and 2) water vapor formed by 
evaporation and transpiration from the land surface of the region. This second 
source is tantamount to the local recycling of precipitation water within the 
region. 

The sources of precipitation water are indicated schematically in 
Figure 2.1, which shows a simplified model of the atmospheric moisture fluxes 
over a land region. The horizontal arrows indicate the advective flux of water 
vapor into and out of the atmospheric control volume; W is the amount of 
water vapor contained in the control volume; E is the net evapotranspiration 
from the underlying land surface; and P is the net precipitation onto the land 
surface. The two sources of precipitation are indicated by the two small 
branches that join to form the larger arrow labeled P; P m is precipitation of 
local (evaporative) origin, and P is that of advective origin. The arrow 
labeled E splits into two branches, indicating that a certain fraction of the 
locally evaporated or transpired water is not returned to the land surface as 
precipitation, but joins the atmospheric vapor reservoir and is advected out of 
the control volume. 

In discussing precipitation recycling, it is important to clarify the term 
"recycling." In this study, recycled precipitation is defined as water that 
evaporates from the land surface within a specified control volume and falls 
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Figure 2.1: Simplified model of the atmospheric moisture fluxes over a land 

region. P m and P a are precipitation of local evaporative and 

advective origin, respectively. 
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again as precipitation within the same control volume. The remainder of the 
precipitation falling inside the control volume is considered to be of advective 
origin, without regard to whether it most recently evaporated from an ocean 
surface or from a land surface outside the control volume. With recycling 
defined in this manner, the distinction between recycled and advected precipi- 
tation is not necessarily the same as the distinction between precipitation 
water of terrestrial origin (i.e., most recently evaporated from the land surface) 
and oceanic origin (i.e., most recently evaporated from the ocean surface). 
Advected precipitation water is equivalent to oceanic precipitation water if and 
only if the inflow boundary of the control volume lies along the coastline. For 
an inland region, it is important to assure compatibility of definitions before 
comparing different estimates of the sources of the region’s precipitation. 

The extent to which a region’s precipitation is supplied by locally- 
evaporated water is an indicator of the importance of land surface processes in 
the water balance of the region and may also be an indicator of general 
climatic sensitivity to land surface change. The recycling process is a 
potentially significant climate feedback mechanism and land surface-atmosphere 
interaction, which may contribute to the persistence and intensification of 
droughts (Rodriguez-Iturbe, Entekhabi and Bras, 1991a). The objective of the 
work described in this Chapter is to obtain quantitative estimates of the 
degree to which land-atmosphere moisture recycling is active over several large 
continental regions. 

2.2 Literature Review 

Until the late 1930’s, according to Benton et al. (1950), it was generally 
believed that the water precipitated over the continents was directly derived 
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from moisture evaporated from the continents. Under this view, the flux of 
water in the atmosphere from the oceans to the continents was assumed to be 
only that required to replenish the loss of water from the continents by runoff. 
Apparently, there was no conception that vast quantities of incoming water 
vapor from the oceans could pass over continental regions without 
precipitating, or that water vapor which evaporated from the continent could 
be carried away by the wind as atmospheric runoff. 

Benton et al. critiqued the theory that precipitation is largely 
land-derived and described developments leading to a more correct description 
of the role of the atmosphere in the hydrologic cycle. They emphasized two 
concepts: First, vertical motions in the atmosphere are necessary in order to 

produce significant precipitation, therefore increasing the water vapor in the air 
over a region does not necessarily increase precipitation. Second, the 
atmosphere is in continuous motion, carrying very large quantities of water 
vapor across the continents from the oceans; moisture added to the atmosphere 
by evaporation and transpiration may travel hundreds or thousands of miles 
before being reprecipitated. 

In the same paper, Benton et al. estimated the relative contributions of 
the sources of precipitation over the Mississippi watershed. They distinguished 
between precipitation events from, and evaporation into, maritime and 
continental air masses, and estimated that not more than ten percent of the 
total precipitation of the Mississippi watershed is moisture having a land 
source within the watershed, while at least 90 percent is external in origin. 

McDonald (1962) argued against what he called "the evaporation-precipi- 
tation fallacy," that is, the idea that local water shortages might be alleviated 
by creating open areas from which water could evaporate and enhance local 
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precipitation. In arguments parallel to those of Benton et al., he attributed 
the fallacy to misconceptions concerning the magnitude of and distance scales 
involved in atmospheric water vapor transport. 

Budyko (1974, pp. 239-243) used a simple one-dimensional model (which 
is described in Section 2.3) to estimate, for the European U.S.S.R., the 
contributions of locally evaporated and advected moisture to regional 
precipitation. Budyko defined /3 as the ratio of total precipitation to 
precipitation due to moisture advected from outside the region. On an annual 
basis, only about 10 percent of the precipitation in the European U.S.S.R. is 
of local origin, according to Budyko’s estimate; on a monthly basis, the 
estimated contribution of locally-evaporated water ranged from four percent in 
October to 18 percent in April and May. Budyko stated, "Even on the most 
extensive continents, where the relative role of local evaporation is the 
greatest, as calculations show, the main portion of precipitation is formed from 
water vapor of external origin, not local." Shiklomanov (1989) used aerological 
data to enhance Budyko’s model, and obtained monthly estimates of the local 
fraction of precipitation in the European U.S.S.R. ranging from .1 percent in 
January and February to 20 percent in June. 

The significance of the oceans as sources of precipitation water to the 
continents is no longer disputed. This understanding of the atmospheric 
branch of the hydrologic cycle has been strengthened by analysis and 
observation during the last half century. However, during the same period, 
studies have shown that land regions can also be significant sources of water 
vapor to the atmosphere. 

Thornthwaite (1946) observed that the average July evapotranspiration in 
the eastern United States (between 5.5 and 6.0 inches) exceeded the July 


52 



evaporation from the water surface in the Gulf of Mexico and the Caribbean 

(about 4.0 inches). Benton and Estoque (1954) showed that the North 

American continent, as a whole, is a source of moisture to the atmosphere in 

summer (May-August) and a sink during the rest of the year. Starr and 

Peixoto (1958) observed strong maxima of vapor flux divergence over three 

separate arid regions, implying a large excess of evaporation over precipitation 

in these areas. Stidd (1967) stated, 

If we consider that in most areas the soil will be drier in the fall 
than in the spring and that the difference will represent a net loss of 
moisture to the atmosphere (in the absence of runoff), it is not hard 
to understand that land areas in general furnish a source of moisture 
to the atmosphere in summer and that the oceans provide a sink, so 
that, in effect, the hydrologic cycle is reversed. 

If the land surface acts as a significant source of moisture to the 
atmosphere during certain seasons, then it is reasonable to expect that the 
balance of locally-evaporated versus advected precipitation might be shifted 
during those seasons (assuming the presence of the dynamical processes 
required to produce precipitation at all). Furthermore, it is reasonable to 
expect that alterations in the storage properties or to the recharge-discharge 
cycle of the land surface and subsurface might affect the precipitation regime. 

Stidd (1967) studied local modifications in climate following large-scale 
irrigation development in the Columbia River Basin. His analysis showed an 
increase in July and August precipitation, extending to several thousand square 
kilometers around the project. Fowler and Helvey (1974) re-examined the 
topic by alternate procedures and concluded that a large increase in 
precipitation due to nearby irrigation was improbable, and that Stidd’s claim 
of a significant increase in precipitation did not appear sound. Stidd replied 
(1975) with statistical analysis supporting his earlier claim. He argued that 
Fowler and Helvey examined the problem on too small a scale, citing the fact 
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that the irrigation project lies near the center of the drainage basin, and that 
the additional moisture added to the air by irrigation would be expected to 
precipitate not immediately adjacent to the irrigation project, but in the 
foothills downwind, as the moist surface air is carried over the surrounding 
mountains. This explanation is in agreement with the statements of Benton et 
al. (1950) and McDonald (1962) that increasing the water vapor content of the 
air will increase precipitation only if an uplift mechanism exists. 

Natural or anthropogenic changes that enhance (or inhibit) convection 
could alter the contribution of local moisture to precipitation. Dettwiller and 
Chagnon (1976) found upward trends in the warm season rainfall at Paris, 
St. Louis and Chicago, and suggested that the urban heat island contributes to 
larger and more intense shower clouds. Schickedanz (1976, pp. 99-100) found 
evidence for irrigation-related increases in summer rainfall in several areas of 
the U.S. Great Plains; he hypothesized that 

. . . any increased rainfall does not come directly from the increased 
atmospheric moisture alone, but by thermodynamic and physical side 
effects produced by the presence of a cool, moist dome over the 
irrigated area. 

Anthes (1983) speculated that if human activities can inadvertently affect 
precipitation, perhaps humans could make intentional changes in surface 
characteristics so as to modify precipitation in a useful way. He proposed 
planting bands of vegetation in semi-arid regions, hypothesizing that the bands 
could increase convective precipitation through three mechanisms: increased 

low-level moist static energy, the generation of mesoscale circulations, and 
increased atmospheric water vapor. 

Lettau et al. (1979) used climatonomic methods (Lettau, 1969) to 
quantify precipitation recycling in the Amazon River Basin, in an effort to 
assess the possible climatic consequences of large-scale deforestation in 
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the region. Among their results were estimates of the ratio of total regional 
precipitation to the regional flushing of precipitable water of direct oceanic 
origin (7). In the six 5' wide regions that they analyzed, 7 increased with 
distance downwind of the coast, reaching a maximum value of 1.884. This 
value of 7 indicates that 47% of the precipitation falling on the subregion at 
75° W has been most recently evaporated from the continent, both upwind of 
and within that particular subregion. It should be noted that Lettau et al.’s 
comparison of 7 with Budyko’s 0 for this inland region is not valid; 0 
compares total precipitation to that of advective origin, where the advective 
term includes moisture evaporated from the land lying between the coast and 
the inland region’s boundary in addition to water of oceanic origin, whereas 7 
compares total precipitation to that of oceanic origin only. 

Salati et al. (1979) studied the inland gradient of the oxygen-18 content 
of precipitation in the Amazon Basin. For yearly— averaged station data, the 
gradient d(<5 18 0)/dx is much smaller than in other continental areas. The 
small inland decrease in isotopic content of precipitation is an indication that 
significant amounts of moisture are added to the air mass as it passes over 
the region, and that this re-evaporated moisture is important in precipitation 
falling on the region; this effect varies with the season and the location. 
Salati et al. did not give a numerical estimate of precipitation recycling, but 
they cited a study by Marques et al. (1972), which found that inflowing 
moisture accounted for only 52 percent of the precipitation between Belem and 
Manaus. 

In atmospheric General Circulation Models (GCMs), water can be tagged 
according to its evaporation site and traced in order to determine the relative 
contributions of different evaporative sources to a region’s precipitation. 
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Joussame et al. (1986) conducted such an experiment for the month of July, 
using the Laboratoire de Meteorologie Dynamique GCM; they produced global 
charts showing the influence of ocean and continental source regions in 
continental precipitation. Koster et al. (1986) used the NASA/Goddard 
Institute for Space Studies GCM for a similar experiment, for all four seasons. 
Koster et al. (1988) also studied the inverse problem, that is, the characteristic 
distances and directions traveled by water evaporating from a source region. 
Shukla and Mintz (1982) conducted a sensitivity experiment with the Goddard 
Laboratory for Atmospheric Sciences GCM, in which two different constraints 
were placed upon the land-surface evapotranspiration: in the first case, no 

evapotranspiration was allowed, and in the second case, evapotranspiration was 
set equal to the model-calculated potential evapotranspiration. They found 
that "land-surface evapotranspiration has a large influence on the precipitation, 
temperature, and motion fields of the atmosphere." 

The importance of land regions in supplying moisture to the atmosphere 
has been well demonstrated. There are indications that land— evaporated 
moisture supplies a large fraction of the precipitation for some continental 
regions. Questions remain concerning the contribution of local evaporation to 
local precipitation. Several quantitative estimates have been made, and there 
is general agreement on the physical factors necessary for such a recycling 
process to occur. 

Recent research by Rodriguez-Iturbe et al. (1990a, 1990b, 1991) and by 
Entekhabi et al. (1991) has shown that the close coupling between the land 
surface and the atmosphere in continental-type climates helps to explain the 
statistical structures exhibited by climatic variables. In their work, the 
coupling is explicitly represented by including the precipitation recycling 
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mechanism in the land surface water balance equation. The resulting 
functional dependence of precipitation on soil moisture, together with the 
stochastic nature of soil and climate parameters, results in a bimodal 
probability distribution of soil moisture states, corresponding to drought and 
pluvial conditions. Further, the dynamics of the soil moisture equation can 
exhibit fixed point, limit cycle and chaotic behavior. 

The importance of precipitation recycling as an index of land 
surface-atmosphere coupling is a compelling motivation for further study of 
this process and for quantification of the degree to which it is active over a 
variety of continental regions. 

2.3 Recycling Model 


Budyko’s Model 

Budyko (1974) considered a land region having length scale i, average 
precipitation P, and average evapotranspiration E. The average precipitation 
is composed of an advective portion (P ) and a local evaporative portion 

3 / 

(P m ), i-e., 


P 


P + P^ 
a m 


( 2 . 1 ) 


The water vapor content of the air moving across the region is also composed 
of an advected and an evaporated portion. As shown in Figure 2.2, air enters 
the region at velocity u and with moisture content w. The vertical flux 
quantities- -P , P and E--are treated as constants equal to their average 
values; therefore the locally-evaporated moisture content of the air increases 
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Local ly-evapor at ed Advected moisture 



Figure 2.2: Fractions of advected and locally-evaporated moisture in the air 

as it moves across a region of length l. 
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linearly, and the advected moisture content decreases linearly, as the air moves 
across the region. It follows that the average horizontal flux of advected 
moisture over the region is 

<p a 

Q a = w » - ( 2 2 ) 

and the average horizontal flux of locally-evaporated moisture is 

«m - & ~ P J • (2-3) 


The atmosphere is assumed to be fully mixed, so that the ratio of advected to 
locally-evaporated water falling as precipitation is equal to the ratio of 
advected to evaporated moisture present in the air. That is 



l P 


wu — 


I< E - p m) 


(2.4) 


Equations (2.4) and (2.1) are solved for the ratio of total to advected 
precipitation, which is Budyko’s recycling coefficient (/?): 

0 = F- = 1 + Us • (2-S) 

a 

The second term on the right-hand side of (2.5) is the ratio of 
locally-evaporated to advected precipitation. In this model the relative 
contribution of locally-evaporated water is directly proportional to the 
evaporation rate and the length scale of the region, and inversely proportional 
to the rate at which external moisture enters the region. The inverse of 0 is 
the fraction of precipitation due to advective origin, i.e., 



( 2 . 6 ) 
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and the fraction of precipitation due to local origin is 



Shiklomanov’s (1989) estimates of the recycling coefficient made use of 
the improved availability of aerological data, rather than estimates of average 
moisture content and wind speed. The equation or equations that Shiklomanov 
used to compute the recycling coefficients are not presented in his paper; 
however, the tabulations indicate that the procedure used was based upon 
extending Budyko’s arguments to a two— dimensional land region, with moisture 
influx and efflux through the sides of an atmospheric control volume. As 
shown in Figure 2.3, the boundary of the region consists of a segment or set 
of segments (7. ) across which the atmospheric moisture flux is inward, and a 
segment or set of segments (7 ) across which the flux is outward. 
Integrating Q • separately over the two parts of the boundary gives F + , 
the influx, and F“, the efflux, of atmospheric moisture through the sides of 
the control volume: 

F + e - | <3 • n^ d 7 (2.8a) 

7- 

'in 

and 

F- i ( 3 • n 7 d7 (2.8b) 


in which n is the outward unit normal vector. 

7 

By definition, F + contains only advective moisture; F 
advected moisture that remains after P 0 is removed, as well as 

cl 


local origin: 


contains the 
moisture of 
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Figure 2.3: Plan view of a region, 

across which the atmosphei 

outward (7 V 

' out' 
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F = ( F+ - P a A) + (E - P m )A 


(2.9) 


in which A is the area of the region, and P a , P m and E are the rates of 
advective and local precipitation and evaporation per unit area. If the average 
horizontal flux of advected moisture over the region is taken as the arithmetic 
mean of the incoming and outgoing advective moisture, then 


F+ + 


^a = 


( F ‘ 


- p a A ) 


, p a A 

= P+ -Hr 


( 2 . 10 ) 


Likewise, the average horizontal flux of locally-evaporated moisture is the 
mean of the incoming and outgoing local portions, 


0 + ( E - p m ) A ( E " p m ) A 
Q m = 2 5 


( 2 . 11 ) 


Invoking, as before, the assumption of a well-mixed atmosphere, the recycling 
ratio is 


Equation (2.12) 


_ = 1 + M 

a 2F 1 


is identical to (2.5), with F 


+ ‘ 

replacing wu, 


( 2 . 12 ) 


and A replacing t 


Application of the Model in the Present Study 

Depending upon the geometry of the region and the spatial variability of 
3, F + as defined in (2.8a) may not be strictly analogous to wu in Budyko’s 
one— dimensional model. Figure 2.4 shows an example of a situation in which 
the analogy is weak; the mean wind blows at an angle, ot, to the axis of the 
region. The moisture influx across segment 7 j is greater than that across 
segment y however, the air entering across 7 2 travels greater distances across 
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Figure 2.4: Plan view of a region in which F + , the net moisture influx, is 

not strictly analogous to wu in Budyko’s model. 
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the region. A more representative value of the moisture influx is an average 
value weighted according to the distance traveled across the region. 
Alternatively, one might delineate the region parallel to the direction of flow; 
however, the latter approach is undesirable for two reasons. First, the 
direction of flow may change with the seasons, whereas the study region 
should remain fixed. Second, the need for interpolation is minimized if regions 
are defined by connected nodes on the latitude-longitude grid of the data set. 
The following procedure is used in the present study to determine effective 
values of i and wu in Budyko’s original equation: 

A simple rectangular region is defined by connected nodes on the grid of 
the data set. A spatial average of 4 (3), is computed by areally averaging 
the components of ^ over the region. A unit direction vector is defined, 

u = (2.13) 

i <a> ii 

which is, essentially, the moisture-weighted average direction of air flow across 
the region. This direction establishes one or two sides of the rectangle as the 
influx boundary, and one or two as the efflux boundary. 

The effective length, t p is the average distance across the region in the 
direction of u. The expression for / is derived in Appendix A; the result is 
as follows: 

l = Ax = A (2.14) 

eff Ay | cos 6 | + Ax|sin0| Ay|cos0| + Ax|sin0| 

where 6 is the angle formed by u and the positive x axis, Ay is the 
south-north length of the region, and Ax is the mean west-east length of the 
region, 

Ax 5 ■ (2-15) 
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The definition of Ax arises from the latitude-dependent relation of degrees 
longitude and west-east distance on the globe. 

The moisture influx term, (wu| ff) is computed as the length-weighted 
average of the moisture flux vector component parallel to u across the influx 
boundary, i.e., 

'7 

C (3 • u) 4 7 ) sin a d 7 

(wu) eff = ^ j (2.16) 


where 7 a and 7 ^ are the endpoints of the influx boundary and a is the angle 
formed by u and cFy (the derivation of this expression is given in 
Appendix A). 

The recycling model, incorporating the effective parameters, is 


fi = 



1 + 



2 (wu) rff ' 


(2.17) 


Equation (2.17) is formally identical to (2.5); however, the spatial variability 
of the aerological data is incorporated through use of the effective length and 
effective moisture influx terms. The recycling coefficient, /?, as estimated by 
(2.17), is representative of the region only to the extent that u is 
representative. 


2.4 Estimation Techniques and Study Regions 

The ratio of local to total precipitation, P M /P, was estimated using (2.7) 
and (2.17). The effective length and effective moisture influx terms were 
computed from the GFDL data as described in Appendix A. 
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(2.18) 


Evapotranspiration was computed from (1.7), neglecting i e > 

(E) % (P) + ^ | 3-n 7 d7 = (P) + 

in which all quantities now refer to mean monthly values, and the overbar is 
dropped for simplicity. The rationale for estimating evapotranspiration as a 
residual is that the atmospheric water balance reflects actual 

evapotranspiration, whereas pan evaporation data estimate potential 
evapotranspiration, and model equations for actual evapotranspiration generally 
require estimates of soil and vegetation parameters that would be difficult to 
obtain for all regions, and on continental scales. 

The divergence term, (V-3), was computed from the GFDL data as 
described in Section 1.4. Regional average precipitation, (P), was computed 
by the Thiessen polygon method from the NCAR station data as follows. 
First, (P) was computed for each month from May 1963 through April 1973, 
then the ten January values, ten February values, and so on, were averaged to 
give a set of monthly average precipitation figures representative of the 
10-year period. Estimation by (2.17) is referred to as the Modified 1-D 
method hereafter. 

The ratio P m /P was estimated using (2.12), in order to make a 
comparison of the two methods and to evaluate the claim that (wu) eff is more 
representative of the moisture influx than F + . Estimation by (2.12) is referred 
to as the 2-D method. The 2-D method is computationally simpler than the 
Modified 1-D method because the terms F + and F can be incrementally 
summed as trapezoidal rule integration proceeds along the boundary, by adding 
each f^ to the appropriate sum according to its sign. The Modified 1-D 
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method is currently developed only for rectangular study regions, which limits 
its usefulness somewhat. 

Besides the rect angularity requirement, there were no constraints on the 
selection of continental regions for analysis. As explained before, the method 
is general, in that the distinction is made between local and advected 
precipitation, rather than terrestrial and oceanic precipitation. The method 
can, therefore, be applied to regions in the interior of continents as well as 
coastal regions. Regions were selected on the basis of hydrological interest and 
the existence of former recycling estimates with which to compare the results. 

The four regions selected for analysis are shown in Figure 2.5. The 
North American region lies east of the Continental Divide and within the 
Mississippi basin (except for the northeast corner, which lies in the Great 
Lakes drainage). Several estimates of local recycling or of the land-evaporated 
fraction of precipitation have been made for this region, which receives 
atmospheric moisture from the west and the south. The European region 
corresponds roughly to the European part of the U.S.S.R., a region for which 
monthly water budgets and recycling ratios have been reported by Budyko 
(1984) and Shiklomanov (1989). The South American region contains most of 
the Amazon river system and rainforest. This region is currently the focus of 
intense hydroclimatological interest, due to the threat of massive deforestation, 
with possible ramifications for the region’s climate, and the climate system of 
the Earth as a whole. The African region lies south of the Sahara Desert and 
includes most of the drainage basin of the Niger River. This region is 
threatened by drought and desertification, and its annual cycle of atmospheric 
moisture supply exhibits unique characteristics. 
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Figure 2.5: Regions selected for precipitation recycling analysis. The analysis 

technique dictated a rectangular shape; otherwise, the regions 
were chosen on the basis of hydrological and historical interest. 
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2.5 Results 


Eurasian Region 

The results for the Eurasian region are presented in graphic form in 

Figures 2.6 through 2.8, and in tabular form in Tables 2.1 and 2.2. By the 
Modified 1-D method, the greatest recycled fraction of precipitation (Figure 

2.6) is 0.39 in June, and the lowest is 0.0 in February. By the 2-D method, 

the largest P m /P is also in June, but somewhat lower than the Modified 1-D 
result. Figure 2.7 shows separate plots of the various terms used in 
calculating P /P. The annual march is as expected for atmospheric moisture 
convergence, precipitation, and evaporation. The fourth plot in Figure 2.7 
compares the influx terms used in the two estimation methods. Equation 

(2.17) can be rewritten as 


0 


1 + 


E 


(wu) 


eff 


' e ff 


(2.19) 


and (2.12) as 


0=1 + 


2 


E 

F + 


( 2 . 20 ) 


The terms in the denominator both represent normalized influx terms and have 
units of length. For the Eurasian region, (wu) eff /^ ff is consistently less than 
F + /A except in July. The difference is greatest in May and June. Figure 2.8 
is included to illustrate the reasons for such differences. Figure 2.8 is a set of 
bar graphs of the boundary flux of water vapor, as introduced in Section 1.4. 
The pattern in May recalls Figure 2.4, the constructed example of a situation 
on which F + is not perfectly analogous to Budyko’s wu. The bars in Fig. 2.8 
do not represent vector components, as do the lines in Fig. 2.4; nonetheless, 
the same arguments apply: (wu) eff /^ eff is less than F + /A because the former 
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Figure 2.6: 


Eurasian Region. The ratio of locally-evaporated to total 
precipitation, as estimated by the Modified 1-D and the 2 L> 
method. 
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Figure 2.7: Eurasian Region. Terms used in estimating P^/P by two 

methods. Atmospheric moisture convergence (-V-3), Precipitation 
(P) and Evaporation (E) are common to both methods; wu II 

efr eff 

is used in the Modified 1-D method, F + /A in the 2-D method. 
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Figure 2.8: Eurasian Region. Bar graphs of the moisture flux across 

segments of the boundary. Each bar represents the integral of 

the normal component of <3 along the segment at the base of the 
bar. 
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Table 2.1 Eurasian Region. Components in the estimation of P /P. 
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Table 2.2 Comparison of estimates of the recycled fraction of precipitation and intermediate terms 
for the Eurasian region. 
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(1) Assuming area = 5.76E+6 km A 2 (2) Assuming 1 — 2300 km 



weights the flux at the west end of the southern boundary more heavily than 
that at the east end, whereas the latter gives equal weight to all the influx. 
In June, the direction of moisture transport across the region is nearly due 
eastward; therefore (wu) eff /^ ff neglects the small influxes over the northern 
and southern boundaries, whereas F + includes them. The computations for 
both the Modified 1-D and 2-D methods are summarized in Table 2.1. 

In Table 2.2, the intermediate calculations and the final P /P results 

m' 

are compared to the work of Budyko (1974) and Shiklomanov (1989). This 
study is in general agreement with their results in all but the summer months, 
lending further support to the claim that external moisture sources dominate 
in supplying precipitation to this region for much of the year. The differences 
in the summer months may simply be due to use of different study regions 
(the region considered in this study does not include the entire European 
U.S.S.R, which was the study region for the former work) and different data 
sets. 

North American Region 

The results for the North American region are presented in Figures 2.9 
through 2.11 and in Table 2.3. The order of presentation is the same as for 
the Eurasian region. Referring to Figure 2.9, P m /P as computed by the 
modified 1-D method lies between 0.2 and 0.4, with the notable exception of 
0.63 in July and August. The minimum P m /P by the Modified 1-D method 
occurs in April. By the 2-D method, P m /P lies between 0.15 and 0.35, with 
minima in April and December; in July and August P m /P by the 2-D method 
is half as large as by the Modified 1-D method. Referring to Figure 2.10, the 
annual march of -V-<3 computed for this region is highly questionable. The 
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Figure 2.9: North American Region. The ratio of locally-evaporated to total 

precipitation, as estimated by the Modified 1-D and the 2-D 
method. 
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Figure 2.10: North American Region. Terms used in estimating P m /P by two 

methods. Atmospheric moisture convergence (-V-cJ), 
Precipitation (P) and Evaporation (E) are common to both 

methods; wu eff /^ eff is used in the Modified 1-D method, F + /A in 

the 2-D method. 
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LOCATION MAP 



Figure 2.11: North American Region. Bar graphs of the moisture flux across 

segments of the boundary. Each bar represents the integral of 

the normal component of ^ along the segment at the base of the 
bar. 
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Table 2.3 North American Region. Components in the estimation of P /P. 
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top frame of Figure 2.10 shows atmospheric moisture convergence into this 

region only during April, and of very low magnitude. If -V-<3 is 

systematically underestimated, then E is overestimated, which would inflate the 

P /P estimates, 
m' 

For this region, F + /A is significantly higher than ( wu ) eff /^ eff f° r most 
the year, which results from the fact that a large part of the moisture influx 
occurs across the southern boundary (Figure 2.11), while at the same time the 
average direction of transport, 6 in Table 2.3, is toward the east or northeast. 
As a result, the magnitude of (wu) eff is only slightly affected by the largest 
boundary influxes. In particular, for the months of July and August, the 
boundaries that are defined by 6 as "influx” boundaries actually contain 
segments where the flux is outward; furthermore, these outward fluxes (in the 
southwest corner), although small, are the most heavily weighted terms in the 
computation of ( wu ) eff - 

For two GCM grid squares roughly corresponding to the region 
considered here, Roster et al.’s (1986) GCM tracer test showed that the 
percentage contributions of the North American source region to local 
precipitation were 36.4 and 41.6 percent in winter, and 80.1 and 61.9 in 
summer, for the western and eastern grid squares, respectively. The 

comparison is not exact, because Roster et al.’s North American source region 
is the entire continent, whereas this study considers only the evaporation and 
transpiration from within the study region. However, it is interesting to note 
the qualitative agreement between the estimates. 

Both estimation methods give P m /P considerably greater than Benton et 
al.’s (1950) estimate that 10 percent or less of precipitation on the Mississippi 
watershed has its evaporative source within the region. Their estimate was 
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based on distinguishing between continental and maritime air masses, with the 
assumption that evapotranspiration is greater into the former, and precipitation 
greater from the latter. The methods in the present study lump all air into 
the same averaging process, essentially ignoring any correlation between 
evapotranspiration or precipitation and air mass source. 

South American Region 

The results of the South American region are presented in Figures 2.12 
through 2.14 and in Table 2.4. The order of presentation is the same as for 
the preceding regions. Referring to Figure 2.12, P m /P by the Modified 1-D 
method is less than 0.3 for all twelve months, with maxima of 0.27 in October 
and November, a minimum of 0.13 in June, and values near 0.15 from July 
through September. The results by the 2-D method are very similar, with 
slightly higher values in the Southern Hemisphere summer, reaching 0.32 in 
December. Referring to Figure 2.13, divergence calculations show the region to 
be a sink of atmospheric moisture throughout the year, but particularly in the 
first four months of the year, corresponding to the southernmost excursion of 
the Intertropical Convergence Zone (ITCZ) during the Southern hemisphere 
(SH) summer. The season of least recycling is also the season of least 
precipitation and least evapotranspiration, the SH winter. Even then, the 
magnitudes of E and P are quite high, as can be seen by comparing Figure 
2.13 to Figure 2.7 and Figure 2.11 (the latter with caution, noting the 
probable overestimation of E for the North American region.) 

The boundary flux pattern (Figure 2.14) illustrates a conceptual problem 
in the Modified 1-D method. Consider the month of July. During July, the 
average direction of moisture flux was determined to be 168° (i.e., from the 
southeast), defining the southern and the eastern boundaries as influx 
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Figure 2.12: South American Region. The ratio of locally-evaporated to total 

precipitation, as estimated by the Modified 1-D and the 2-D 
method. 
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Figure 2.13: South American Region. Terms used in estimating P m /P by 

two methods. Atmospheric moisture convergence 
Precipitation (P) and Evaporation (E) are common to both 

methods; wu tff /^ ff is used in the Modified 1-D method, F + /A in 

the 2-D method. 


83 






LOCATION MAP 



-75 -5C 

Longitude, Degrees 



Figure 2.14: South American Region. Bar graphs of the moisture flux across 

segments of the boundary. Each bar represents the integral of 

the normal component of ^ along the segment at the base of the 
bar. 
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Table 2.4 South American Region. Components in the estimation of P / P. 
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53.5 212.9 159.4 -168.5 2.18E+06 4.86E+08 222.8 1.36 0.74 0.26 168.7 1.47 0.68 0.32 



boundaries. However, during July there is a bifurcation of the mean mass 
flow over the region, with the result that the southern boundary is actually an 
efflux boundary, and 3 across the boundary contributes negatively to (wu) ef{ 
(Eq. 2.16). Nonetheless, the months when the bifurcation is strongest are also 
the months when F"*" and (wu)^ are most similar, indicating the need for 
further analysis and comparison of the estimation methods. 

Thirty percent recycling for the region as a whole is compatible with 

Lettau et al.’s (1979) estimate that 47% of the precipitation falling at 75° W is 

of land-evaporated origin. If the contribution of local moisture varies linearly 

(Figure 2.2), then the fraction of locally evaporated water at x = l will be 

greater than the average value. This study indicates that, on average, more 

than 70 percent of the region’s precipitation is provided by inflowing moisture. 

This estimate is considerably higher than that of Marques et al. (1972), that 

inflowing moisture accounts for 52% of the precipitation between Belem and 

Manaus. Belem, at the mouths of the Amazon, lies to the east of the present 

study region, and Manaus (longitude 60° W), east of the center of the study 

region. By the same argument as above, the Modified 1-D model predicts 

that the advective fraction of precipitation for x less than 1 / 2 will be greater 

than the average value; that is, P a /P for the longitudes between Belem and 

Manaus is between 70 and 100 percent, or in terms of the local contribution, 

P / P is between 0 and 30 percent. Subdividing the study region would 
nr 

provide a more adequate basis for comparison. 

African Region 

The results for the African region are presented in Figures 2.15 through 
2.17 and in Table 2.5. The order of presentation is the same as for the 

preceding regions. Referring to Figure 2.15, two peaks appear in the annual 
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Figure 2.15: African Region. The ratio of locally-evaporated to total 

precipitation, as estimated by the Modified 1-D and the 2-D 
method. 
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Figure 2.16: African Region. Terms used in estimating P m /P by two 

methods. Atmospheric moisture convergence (-V -<3), 

Precipitation (P) and Evaporation (E) are common to both 

methods; wu „/l ,, is used in the Modified 1— D method, F + / A in 
’ err eff 

the 2-D method. 
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LOCATION MAP 



OCT NOV DEC 

Figure 2.17: African Region. Bar graphs of the moisture flux across segments 

of the boundary. Each bar represents the integral of the normal 

component of <3 along the segment at the base of the bar. 
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Table 2.5 African Region. Components in the estimation of P m /P. 
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march of P m /P: 0.52 in February and 0.64 in August by the Modified 1-D 

method; and 0.41 in February-March and 0.48 in August by the 2-D method. 
The February-March peak corresponds to fairly high E and low P in those 
months, while the July-August peak corresponds to very high E and P (Figure 
2.16). Analogous to the South American region, large differences in the two 
influx terms occur in months of strong convergence; the month of June is the 
best example (Figure 2.17). The convergent boundary flux during June 

corresponds to the northernmost excursion of the ITCZ during the northern 
hemisphere summer. During June, according to calculations, the average 

direction of moisture flux is 178°, nearly due westward; as a result (wu) ff 
virtually ignores the influx over the southern boundary, whereas F + includes 
that influx, as well as the influx across the northern boundary. 

Besides the June convergence noted above, Figure 2.17 illustrates several 
interesting characteristics of the atmospheric water balance for this study 
region. In February and March, the striking difference between the 
inward-directed bars on the northern boundary and the outward directed bars 
on the southern boundary probably reflects the evaporating power of the dry 
air that moves off the Sahara desert and gains moisture as it traverses the 
study region. 

The estimated monthly rates of moisture convergence (Figure 2.16) 
indicate that, except during May and June, the territory defined as the 
African region is a net source of moisture to the atmosphere. The defined 
region includes the expansive lower Niger marshlands, where a large volume of 
the water drained from the Guinea highlands and the tropical upper Niger 
evaporates into the semi-arid surroundings. There is considerable excess of 
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evaporation over precipitation in this region during the months when the 
highland river discharge spreads over the expansive marshland and evaporates 
into the atmosphere. During May and June when the River experiences its 
lowest monthly flow rates (and even dries completely in some years) the 

marshlands disappear and the regional evaporation is accordingly reduced. 
During this season, the precipitation water is mostly derived from advected 
moisture, as evident in Figure 2.15. Furthermore, there will be a net regional 

convergence of moisture > 0) during these months when the sources 

of water vapor lie outside of the defined territory (Entekhabi, personal 
communication, 1991). 

The significant contribution to precipitation of locally evaporated 
moisture in July and August can readily be inferred from Fig. 2.16. 

Precipitation is at its maximum during these months, yet the influx of 
atmospheric water vapor (by either method of estimation) is much less than in 
June and September; the larger part of precipitation must be supplied by local 
evapotranspiration . 

2.6 Discussion 

Both of the precipitation recycling estimation methods used in this study 
are modifications of Budyko’s linear model of the water vapor content of the 
air moving over a land region. The critical assumption of the model is that 
locally-evaporated and advected water vapor molecules condense and fall as 
precipitation proportionally to their respective fraction of the total number of 
vapor molecules contained in the entire atmospheric column. Stidd (1967) 
argued against the fully-mixed assumption: "Since precipitation is normally 

associated with rising currents of air, it is logical to suppose that most of the 
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excess moisture associated with a storm has come recently from a layer of air 
close to the ground." The larger fraction of locally-evaporated moisture in the 
lower layers then would imply a larger fraction of local moisture in convective 
precipitation. On the other hand, if non— precipitating convection mixes the 
troposphere vertically (Paluch, 1979; Taylor and Baker, 1991), the assumption 
of a well-mixed atmospheric moisture reservoir is strengthened for regions and 
seasons in which convection is active. 

A more serious concern, for the present study, is the problem of air mass 

divergence in the GFDL analyzed data. In this study, the North American 

region (where radiosonde data available for assimilation were rather densely 

distributed, when compared with the rest of the globe) appears to be a net 

source of moisture to the atmosphere on an annual basis. The Mississippi 

River and many of its tributaries provide an inward flux across the horizontal 

boundaries of the region; furthermore, the region is underlain by several major 

groundwater aquifers which are known to be in a depletion trend due to 

pumping for agriculture. The former flux is reflected as a positive R. , and 

in’ 

or 

the latter as a negative ^ in (1.10). However, it is doubtful that these 

additional sources of water are reflected in the atmospheric water balance, and 
such speculation is moot while air mass divergence errors distort the accuracy 
of (V-3). In Savijarvi’s (1988) modified GFDL data (Figure 1.2), the 
year-round divergent area fraction of North America appears to be 
substantially reduced from the original GFDL data (Figure 1.1). 

A priority for further work along these lines is to correct the air mass 
divergences. Alestalo’s method would provide corrected values of V-C$ at grid 
nodes. For the present study’s Modified 1-D recycling estimation method, a 
minor adjustment would be required in order to estimate (V-<3) by an area 
integral rather than a boundary flux integral. However, the 2-D method 
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would require corrected values of the components of in order to distinguish 
'y from 7 , and to compute F + . Savijavi’s correction technique provides 

'in out’ 

corrected components of the mean wind, and is more physically based than 
Alestalo’s technique. 

Another assumption in the recycling estimation methods is the use of 

areally-averaged precipitation over arbitrarily-defined regions. Precipitation is 

an extremely variable process. Although time averaging smooths the series at 

a point in space, substantial spatial variation may exist, depending upon the 

topography (and other surface properties) and atmospheric circulation 

characteristics of a region. As an extreme example, the model is invalid if all 

of a region’s precipitation falls near the upwind boundary, before any 

locally-evaporated moisture has been added to the air. In such a situation, 

P / P would be overestimated, because the model assumes that the same 
m' 

amount of total precipitation falls from the air while the fraction of 
locally-contributed moisture in the air increases linearly. 

A claim was made in Section 2.3 that the distance-weighted average 
moisture influx, (wu) (ff , is more representative of the rate of water vapor 
advection into the region than is the total moisture influx F + . The results 
show that the two terms, when appropriately normalized, are generally 
compatible, with certain exceptions. In three regions (all but South America), 
F + /A is consistently greater than ( wu ) c ff/^ e ff South America, the two are 
nearly equal, with the greatest differences occurring in the months when the 
average moisture flux over the region, (^) , is most southward (as indicated by 
0 in Table 2.4). The supposed superiority of (wu) eff is not clearly confirmed. 
The process of resolving C$ parallel to u at the boundaries essentially neglects 
any curvature in the streamlines of 3 over the region, and eliminates a 
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possibly significant portion of the influx, as illustrated in Figure 2.18. 
Another complication arises, as mentioned before, when the sides of the 
rectangular boundary that are designated as inflow segments by the direction 
of u are not truly inflow segments along their entire length (as for the African 
region in December, Figure 2.17 and Table 2.5). The use of (wu) eff introduces 
unforeseen problems, which call into question the general applicability of the 
Modified 1-D method. The simpler 2-D method may be more appropriate to 
the geometry of Q, but more analysis is needed. In the meantime, it is 
probably safer to believe the lower estimate of P m /P in those cases in which 
a large difference exists between the two estimates, as for the North American 
region in July and August, and the African region in January, February, July 
and August. 

The local precipitation ratio, as estimated herein, is diagnostic. The 
values are indicators of the importance of land-atmosphere interaction to 
regional climate, for climatic conditions during the years 1963-1973. As such, 
the estimates are not predictive; for example, it does not follow that reducing 
evapotranspiration to zero in the Eurasian region during June (even if that 
were possible) would reduce June precipitation by 35 percent. Many complex, 
interrelated factors, both internal and external to the region, control the 
region’s precipitation. 

It is clear that a great deal of variety exists in the importance of the 
precipitation recycling process to different regions and at different seasons. 
Therefore, Budyko’s (1974) conclusion that the relative role of local recycling 
is limited, and Lettau et al.’s (1979) conclusion that it is quite important, are 
not contradictory. The reason is that the process depends, not only on the 
region’s length scale and evapotranspiration rate, but also on the presence of 
the physical mechanisms that cause precipitation to occur at all. 
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Figure 2.18: a) Plan view of a region with convergent moisture flux, and 

with constant |3| along the inflow boundaries. By 

symmetry, the average direction of moisture flux (u) is due 
eastward. 

b) The difference between the total boundary integral of influx 
(F + ) and the integral of boundary influx resolved parallel to 
U [(wu) eff Ay]. 
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If a sufficiently high orographic barrier is present at the downwind 
boundary of a region, most of the moisture present in the air, regardless of its 
evaporative source, precipitates on the windward side. By the definition of 
recycling used in this study, the moisture content due to evapotranspiration 
within the region — which lies simultaneously upwind and downhill — falls back 
into the region as recycled precipitation. This mechanism illustrates the 
sensitivity of the recycling notion to the choice of control volume, and sheds 
some light on why increased irrigation could increase rainfall in the Columbia 
River Basin (Stidd, 1967) whereas building a large lake on the southern border 
of Arizona would not increase rainfall in that state (McDonald, 1962). 

The causal connection between evapotranspiration and convective 
precipitation merits further investigation. Adding moisture to the air at low 
levels tends to reduce the stability of the atmosphere by building up a 
reservoir of latent heat in the lower layers (Stidd, 1967). By this 
contribution, increased evapotranspiration could enhance convective 
precipitation. If so, advective precipitation (P ) would likely be increased as 
well as P m , whereas with the orographic mechanism, increased 
evapotranspiration would only increase P m 

In this study, the partition of precipitation between local and advective 
sources has been emphasized rather than the oceanic/terrestrial distinction. 
For studies of the interactions between the land surface and regional climate, 
the internal vs. external supply of water vapor is of interest, whether that 
external supply comes from an ocean or a land surface. As noted before, the 
definitions merge if the control volume is taken as an entire continent. A 
pragmatic reason for the advective vs. local distinction in the present study is 
that the analysis technique currently dictates the use of rectangular regions, 
incompatible with the shapes of the Earth’s continents. 
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The regions examined in this study were somewhat arbitrarily selected. 
They vary widely in size and hydroclimatological characteristics. More 
systematic selection and study of regions is certainly called for, in tandem 
with improvements to the analysis technique to allow the delineation of 
boundaries along natural drainage divides and coastlines. 

Future work should make use of corrections to the GFDL data, and of 
additional data sets. Four-dimensional data assimilation techniques, now 
routinely used in weather forecasting, provide gridded fields of the variables 
required for recycling estimates. The improved temporal resolution provided 
by these analyses (four times daily, at present) would better capture diurnal 
variations (noted by Rasmusson, 1967) and help correct bias in the <3 
computations now based on once-daily/twice-daily observations. In addition, 
the estimation of (P) should be improved, and methods devised to use other 
water balance quantities such as runoff to provide closure and a check on the 
evapotranspiration term. 
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Chapter 3 

Parameterization of a Simple Climate Model 


3.1 Introduction 

Atmospheric General Circulation Models (GCMs) are the hydrologist’s or 
climatologist’s experimental apparatus for large-scale investigations of the 
atmospheric branch of the hydrologic cycle and land surface-atmosphere 
interactions (Henderson-Sellers and McGuffie, 1987). However, certain 
components of GCMs require improvement before the models can be fully 
useful operationally, for regional studies, or for reliable predictions (Eagleson, 
1982; Entekhabi and Eagleson, 1989). One such component is the 

representation of land surface hydrology. 

General circulation models are based on numerical solutions to the basic 
equations of mass, momentum and energy in the atmosphere. The equations 
are solved at nodes of a three-dimensional mesh, with horizontal grid spacing 
on the order of 10 4 to 10 5 km 2 . Any physical process that occurs on scales 

smaller than the grid resolution must be parameterized. Parameterization 

refers to the process of developing equations that capture the integrated effects 
of the smaller scale processes on the scale resolved by the model. Usually the 
equations require numerical values, or parameters, that characterize the process 
at a particular location or time. 

Several problems arise when new parameterizations are developed for 
GCMs. First, repeated simulations are required in order to test the effects of 
the parameterization; the cost of such extensive computation may be 
prohibitive. Second, feedbacks in the GCM climate make it difficult to isolate 
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specific climatic processes and the effects of the candidate parameterization. 

A one-dimensional screening model has been developed to test possible 
modifications of the GCM land surface hydrology parameterization before the 
changes are implemented in the full GCM (Koster et al., 1988; Entekhabi and 
Eagleson, 1989; Koster and Eagleson, 1990). The screening model is a simple 
representation of the interaction among a land surface, an ocean surface, and 
their overlying atmosphere. Because the screening model does not simulate 
three-dimensional dynamics, a simple linear transfer scheme is used to parame- 
terize the lateral convergence of heat and moisture from the atmospheric 
column overlying the ocean to the column overlying the land surface. 

The objectives of the work described in this Chapter are 1) to determine 
whether the simple linear parameterization of heat and moisture convergence 
produces a realistic environment for the testing and screening of candidate land 
surface hydrology parameterizations; and 2) to determine a range of 
appropriate numerical values for the transport parameter required in the 
screening model. The present study focuses on the moisture convergence. The 
vertical distributions of heat and moisture in the atmosphere are similar, and 
advection is the dominant mechanism of lateral transport for both. Therefore, 
the results for moisture should be applicable to heat as well. 

3.2 The Model 

The one-dimensional screening model is described by Koster and Eagleson 
(1990). A brief summary is given here. 

The hydrologic cycle of the model consists of the interaction between two 
dimensionless regions, one land and one ocean. Basically, the land and ocean 
regions are reservoirs with the appropriate characteristics for storage and 
vertical fluxes of heat and water. To each region there corresponds an 
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atmospheric control volume, defining a box with unit horizontal area whose 
base is the land or ocean surface. 

In the long term, the net transport of water vapor from the ocean box 
into the land box must balance the net liquid runoff from the land box to the 
ocean box. A fixed value of the atmospheric convergence will thus dictate the 
climatic response of the surface. Parameterization using a fixed convergence 
would therefore severely distort the land surface’s response due to changes in 
hydrologic parameterization (soil and vegetation type, etc.) and would be 
counter to the purpose of the model, which is to determine the sensitivity of 
the model climate to such changes. 

Therefore, the lateral convergence rates of heat and moisture must be 
represented in such a way that they may adjust to the changing atmospheric 
profiles that are brought about by changes in model specifications or 
parameterizations. A first-order appproximation for the exchange rate is to 

assume it proportional to the difference in concentration of the 

quantity — specific heat or specific humidity — in neighboring boxes. Thus, 
moisture moves from the oceanic to the continental box at a rate specified as 
follows: 




(3.1) 


in which q is the specific humidity at pressure level i; the subscript c or o 
refers to continent or ocean, respectively; and x is the proportionality constant, 
having units of time. A system such as (3.1), in which the outflow is 
proportional to the storage is known in engineering as a linear reservoir. 
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In the screening model, the land and ocean air columns form a loop, in 
that the inflow to one is the outflow from the other. The moisture exchange 
between boxes is followed by evaluation of the equations controlling the 
various vertical fluxes and readjustments of the atmospheric profile in each 
box. 

3.3 Estimation Technique 

The objectives of the study require a comparison between the moisture 
convergence behavior of the model and of real Earth regions. The procedure 
used to make this comparison is now described. 

The transport parameter x is prescribed for the model and controls the 
rate of moisture exchange between the oceanic and continental boxes. 
Different values of x produce different values of moisture convergence onto the 
continent. A series of ordered pairs (x, -V-C$) can be obtained by assigning 
values of x and integrating the model for a sufficient number of time steps to 
obtain a climatic value of -V-3 (where -V-<3 is the net moisture convergence 
integrated through the depth of the atmospheric column, as defined in Chapter 
One). 

In order to evaluate the parameter x f° r continents and oceans on the 
Earth, Equation (3.1) must be modified slightly. In general, geographic 
regions on the Earth differ from the model in several respects. First, there is 
no Earth region where flow cycles between a single continent and a single 
ocean source region; that is, the real-world ocean— to— continent moisture 
exchange is not a closed loop, as it is in the model. Second, the surface areas 
of the geographic ocean and continental regions may or may not be equal. 
Third, the model does not include surface topography whereas in the 
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geographic continental region elevation changes produce differences in surface 
pressure and consequently the mass of the column of air overlying points on 
the surface is not constant over the land region. 

The model equation for moisture convergence is modified to account for 
the differences in geometry by taking the net moisture influx to the 
continental region as proportional to the difference in water vapor loads 
averaged over the regions. The constant of proportionality is x (where the 
symbol ‘ will denote data estimates of *)- For the data, then, 


d <q c > 

“IT - 


<%> - <q c > 

X 


(3.2) 


where the angled brackets indicate an areal average. The average moisture 
content of the air overlying the continental region is 



(3.3) 


in which M c is the total mass of water present in the atmospheric control 
volume, i.e. , 

M c= II ( S ^ iA < 3 '4> 

continent 


and A £ is the total mass of air overlying the continent, i.e., 

A c = JJ j'*dA. (3.5) 

.. . o 

continent 

Similarly, for the ocean region, the total mass of atmospheric water vapor is 
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(3.6) 


M. 


■ II I 


ocean 


& dA 
g 


and the total mass of air overlying the ocean is 


A„ = 


r P S 

f ^dA 


(3.7) 


ocean 


0 


Assuming that the total mass of air remains constant (this is tantamount to 
neglecting changes in surface pressure), Equation (3.2) may be rewritten as 


d M 1 dM c 
= 7T TT 


M M 
o c 

XT ■ 


(3.8) 


or 


dM c _ 7T 0 *“o 
~ 3T - 


M - M. 


(3.9) 


The left-hand side of (3.9) is the net lateral atmospheric inflow from the 
ocean control volume into the land control volume, which corresponds to — 
Thus x may be estimated from the data as follows: 

A. 


X = 


K 0 


M - M. 


(3.10) 


(-v-3) 

The quantities defined in (3.4), (3.5) (3.6) and (3.7) were evaluated from the 
GFDL gridded data using trapezoidal rule integration, and allowing for the 
spherical surface geometry in the area integrals by the method described in 
Appendix A. Equation (3.10) was used to obtain a value of x for each month 
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of the year. The result is a series of twelve ordered pairs £) for the 

observational data. The model (-V-3, x) and the data (-V-3, x) can then be 
compared. It was assumed that intra-annual changes in the dynamics of the 
ocean-continent interaction would be reflected as different values of x for the 
different months. 

3.4 Study Region 

The one-dimensional screening model is an abstraction, by definition. It 
is not meant to represent any particular location on the Earth. Rather, it is 
a simple one-loop hydrologic cycle, operating between two reservoirs, where the 
conditions of the reservoirs and the rate of exchange between them are allowed 
to vary. 

Few continental regions behave like the land air column in the model. 
First, the model continent is a constant sink of atmospheric water whereas, as 
discussed in the preceding chapters, actual land regions may act as moisture 
sources to the atmosphere at certain times of the year. Second, because of 
the three-dimensional, time-varying circulation of the Earth’s atmosphere, it is 
difficult to define a unique ocean source region corresponding to a given 
continental region. Third, the model has no topography, whereas in real 
regions, topography affects many aspects of climate. 

The region selected for estimation of x was equatorial South America. 
This land region exhibits atmospheric moisture convergence throughout the 
year. Furthermore, as shown by the streamlines of moisture flux (Figure 3.1), 
there is a reasonably well-defined ocean source region. Two land regions and 
corresponding ocean regions were delineated; they are shown in Figures 3.2 and 
3.3. The analysis described in Section 3.3 was applied to each continent-ocean 
pair of study regions separately. 
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Figure 3.2: The continent and ocean regions used to estimate the transport 

parameter, from observations: Experiment 1. 
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3.5 Results 

Series of pairs (-V>3, x) for six different runs of the one-dimensional 
screening model were generated (Entekhabi, personal communication, 1991). 
The six runs, "Run a" through "Run f", differed in the settings of several 
fixed parameters of the screening model. The various runs were conducted 
solely to give a range of model environments from which to collect the 
(-V*Q, x) pairs, and did not constitute an attempt to match the model 
climate statistics to the geographic test regions. 

The model settings used in the six runs are listed in Table 3.1; the 
reader is referred to Entekhabi and Eagleson (1989) for descriptions of the 
parameters and their roles in the model. The model was run for perpetual 
July 1 at 10° North latitude, which assured that the solar radiation received 
at the top of the column remained constant, except for the diurnal cycle; 
therefore any changes to the radiation regime that affected the model 
hydrologic cycle were brought about by the model climate and the composition 
of the atmospheric profile. 

The model was integrated for 250 days, for each of 25 prescribed values 
of x» fo r each of the six runs. The values of x ranged from 0.1 to 20 days. 
The runs produced six (-V-3, x) sequences; they are plotted in Figure 3.4 and 
tabulated in Table 3.2. 

The (-V-<3, x) pairs obtained by evaluating (3.10) for the two geographic 
continent-ocean pairs (Testl-Twaterl and Test2-Twater2) are plotted in 
Figure 3.5. The values of (-V-3, x) from the geographic study regions are 
tabulated, together with some intermediate calculations, in Tables 3.3 and 3.4. 
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Table 3.1 


One-dimensional screening model parameters used in Runs a 

through f. [Source: Entekhabi, personal communication, 1991] 


Latitude: 

Simulation Period: 

Initial Soil-Water: 

Cloud Temperature 
Standard Deviation a; 

Land Surface Wind: 

Soil Depth and Textture: 

Vegetation: 

Surface Albedos: 


COMMON PARAMETERS 
10 North 

250 days (75 days spin-up) Perpetual July 1st. 

Saturated 

3 degrees Celsius 
2 m/s 

1.5 meters (5 layers with ticknesses=0. 10, 0.15, 0.35, 0.40, 0.50 
meters). Composition: 319b Sand , 99b Silt , 609b Clay 
[From Zobler and Cory (1984) data-set for the Amazon] 

Full cover with -15 bars wilting point 

259b Ocean, 359b Land 


SIMULATION PARAMETERS 


Run 

Name 

Soil Moisture 
Coeff.Var.f ] 

PPT Fractional 
Wetting [ ] 

Fixed Ocean 
Temp. [ C] 

Ocean Surface 
Wind [m/s] 

runa 

1.0 

0.6 

28 

4 

runb 

0.5 

0.6 

28 

4 

rune 

1.0 

0.6 

30 

4 

rund 

1.0 

0.3 

29 

4 

rune 

1.0 

0.3 

28 

8 

runf 

0.5 

0.3 

28 

8 
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Figure 3.4: x and for six different runs of the screening model. For 

each run, all other model parameters were held fixed while x was 

varied; the model was integrated over 250 days for each (— x) 
pair. 
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Table 3.2 


CHI 


0 .100 

0 . 150 
0.200 
0.300 
0.400 
0 . 500 
0 . 600 
0 . 700 
0 . 800 
0.900 

1 .000 
1.250 
1 . 500 
1 . 750 

2 .000 
2 .500 

3.000 

3 . SCO 

4 . 000 

5 . 000 

6.000 

7 . 000 

8 . 000 

10 . 000 
20 . 000 


Moisture convergence values for screening model runs a through f. 

[Source: Entekhabi, personal communication, 1991.] 


-Div Q 


Run_A 

Run_B 

Run_C 

0. 
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0. 

596 

0 . 

227 

0 . 

534 

0. 

649 

0. 
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1 . 

155 

0. 

901 

1 . 

330 

0 . 
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1 . 
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2 . 

033 

0. 
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0. 
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1 . 
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0. 
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0. 
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0. 
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0. 
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0. 
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0. 
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0. 
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0. 
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0. 
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0. 

.480 
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.245 
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0 , 
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1 . 

. 000 

0 , 
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0, 

.251 

1 . 

. 010 

0 
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0 
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0 , 
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0, 

.257 

0, 

.198 

0 , 

.865 

0, 
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0 

.220 

0 

.883 

0 

.203 

0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 
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0 

.216 

-0 
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0 
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0 

.237 

-0 
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0 
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0 
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0 

.008 

0 
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0 
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0 
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0 
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0 
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0 
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Run D 

Run E 

Run_F 

0.581 

0.948 
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0.813 

1 . 477 
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0.928 
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1 . 921 
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Figure 3.5: Estimates of x and -V • ^ from the observational data for each of 

twelve months. The observations that do not satisfy model 
assumptions are labeled by month. 





Table 3.3 Estimation of x from observational data, Experiment 1. 
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Table 3.4 Estimation of x from observational data, Experiment 2. 
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For the geographic regions, there were several months that did not match the 
construct of the simple one— dimensional model. During November, December, 
and January in Experiment 1, the computations showed negative values of 

, with moisture convergence onto the continent, giving negative 

values of x- During July and August in Experiment 2, the computations 

showed a net divergence of moisture from the continent, despite a positive 
moisture gradient between ocean and continent. This caused x to be negative 
for those months. In the screening model, x is always positive. Therefore, 
for purposes of comparison, only the positive x’® from the data were used. 
The outliers are included and labeled in Figure 3.5 for the sake of 
completeness 

To improve legibility, the observational (-V*^, x) pairs, excluding the 
outliers, are plotted and labeled by month on semilog plots in Figure 3.6. A 
seasonal cycle in x is observed, with lower values occurring in January, 
February, March, and April and highest values in June, July, August, and 
September. The values of x range over four orders of magnitude, from .2 to 
250 days, although most lie between 10 _1 and 10 1 days. The lower values of 
-V-($ correspond roughly to maximum x> while the highest -V*<3 is associated 
with x around 0.5 day, for both regions. 

Finally, the series of (-V-3, x) all six model runs and the 
observational (-V-3, x) for both geographic study regions are compared in 
Figures 3.7 through 3.12. Each figure shows one model run as a solid line, 
plotted with the observational data from Experiment 1 (indicated by filled 
circles) in the top frame and Experiment 2 (indicated by filled triangles) in 
the bottom frame. 


AT o c 
o 


116 


Chi | day | Chi |day| 



-DivQ [mm/day] 


Figure 3.6: Estimates of x and -V*C$ from the observational data, for selected 

months. The observations that did not satisfy model assumptions 
have been excluded. 
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The model runs all exhibit a similar pattern, with the highest values of 
X corresponding to the lowest values of — and a maximum value of — 
associated with a critical value of x; the critical value lies between 0.2 and 1.0 
day for all six runs. In addition, the model runs show a plateau of nearly 
constant low values of — for x greater than a value lying between 2.0 and 
5.0. The Test2 data show a similar peak and plateau structure, whereas the 

Testl data show a peak, but no plateau. 

In the model, with a large prescribed value of x> the oceanic and the 
continental air columns are virutally isolated from one another. In this case, 
there is little mixing of air between the two halves of the model and the 
magnitude of V - (^ will thus remain low. At lower values of prescribed x m 
the model, the mixing is more efficient and the magnitude of will be 
larger, due to greater moisture exchange between the ocean and the land air 
columns. At even lower values of prescribed x in the model, the mixing of 
the two air columns is essentially complete. In this case, the constraint on 
the magnitude of is the moist convection scheme, which vertically 

distributes the heat and mosture input at the surface interface of the air 
column. At very low values of prescribed x> the rates (as a consequence 
of limits on precipitation) will be somewhat reduced and tend toward another 
plateau. 

The model and data curves in Figures 3.7 through 3.12 do not match 
exactly. They were not expected to match, because the plots are not a 
comparison of theoretical prediction to experimental results. The model runs 
were not an attempt to model the Amazon basin. This is a comparison 
between a simple numerical model of the hydrologic cycle in the abstract, and 
a continental region that exhibits one of the most self-contained hydrologic 
cycles observed on the Earth. 
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Concerning the objectives outlined in Section 3.1, two results emerge: 1) 

the response of moisture convergence to changing x in the model follows a 
pattern that is qualitatively the same as the relationship between and x 

for the South American region. In the sense of qualitative agreement the 
convergence parameterization is concluded to be a reasonable representation of 
the complex dynamics of moisture convergence. 2) Values of x for a real 

Earth region lie in the range 10 _1 to 10 3 days, with values above 
approximately 10 days corresponding to very low rates of moisture convergence 
and values of about .5 day corresponding to maximal moisture convergence. 
These results give an appropriate scale for the selection of x for the screening 
model. 


3.6 Discussion 

The two sets of anomalous observational data that were excluded from 
Figures 3.6 through 3.12 are now discussed. In the case of Experiment 1, x 
was negative for the months of December, January, and February. As shown 
in Table 3.3 and Figure 3.13, atmospheric moisture converged onto the 
continental region during these months, despite the fact that the scaled water 


content of the oceanic control volume 



was less than the water 


content of the continent control volume (M c ). In absolute terms, the water 
contents of the two reservoirs were most nearly equal during these months. 
However, the oceanic air mass is greater, due to the combined effect of the 
ocean region’s greater surface area and the topography of the land region. 
The areal mean specific humidity of the selected ocean region was less than 
that of the land region throughout the year; this result may be attributable to 
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Figure 3.13: Region Testl. Bar graphs of the moisture flux across segments 

of the boundary. Each bar represents the integral of the normal 

component of ^ along the segment at the base of the bar. 


126 



numerical error and to data sparsity for the Southern Hemisphere in general, 
and the oceans in particular. 

For Experiment 2, the months of July and August exhibited net 

divergence of moisture from the continental region, rather than convergence, 

A 

despite the fact that M was greater than M for these months. In 

o c 

Figure 3.1, the JJA map shows a strong northward component to the moisture 
transport off the ocean study region Twater2. This northward component is 
reflected in Figure 3.14 by the large outward bars on the northeast boundary 
of the continental region. Southward from the Equator, the boundary flux of 
moisture from the ocean source region decreases linearly, and reverses direction 
at the southeast corner. Because of the direction of net moisture transport in 
this season, the direct input of moisture from the assumed ocean source region 
to the continental region as a whole is not great. The same flow pattern 
applies to Experiment 1 in these months; however, Region Testl (Figure 3.13) 
is truncated above the latitude where the linear decrease in landward moisture 
flux begins. 

It is clear from Figures 3.13 and 3.14 that the continental regions Testl 
and Test2 do not interact uniquely with the Atlantic Ocean source regions 
Twaterl and Twater2. Atmospheric runoff occurs over the western boundary 
into the Pacific Ocean, over the southern boundary onto an adjoining 
continental region, and over the northern boundary into the Caribbean Sea. 
The region is not a simple one inlet-one outlet hydrologic cycle, as was 
assumed when it was chosen to give real-world values of the model parameter 
X - A more appropriate choice might be to include the Caribbean Sea and 
Eastern Pacific as part of the ocean study region in both experiments; 
although these ocean regions generally do not contribute moisture to the 
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Figure 3.14: 


Region Test2. Bar graphs of the moisture flux across segments 
of the boundary. Each bar represents the integral of the normal 

component of ^ along the segment at the base of the bar. 
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continent, they receive efflux from the continent. In the model, the same 
ocean box serves both functions. It is impossible to find such a situation on 
the Earth; however, including a ring of ocean around a continent might come 
closer than isolating a particular piece of ocean. 

It has been demonstrated that the linear reservoir parameterization of 
ocean to continent moisture transport in the screening model produces a 
qualitatively realistic relationship between the transport and the climatic 
moisture convergence, when compared with one of the most straightforward 
continental hydrologic cycles on the Earth. 


129 



REFERENCES 


Alestalo, M., 1983. The atmospheric water vapor budget over Europe. In 
Variations in the Global Water Budget, ed. A. Street-Perrot et al., 
pp. 67-79. Dordrecht: D. Reidel Publ. Co. 

Anthes, R. A., 1984. Enhancement of convective precipitation by mesoscale 
variations in vegetative covering in semiarid regions. Journal of Climate and 
Applied Meteorology, 23: 541-554. 

Benton, G. S., R. T. Blackburn and V. 0. Snead, 1950. The role of the 
atmosphere in the hydrologic cycle. Transactions, American Geophysical 
Union, 31: 61-73. 

Benton, G. S. and M. A. Estoque, 1954. Water vapor transfer over the North 
American continent. Journal of Meteorology, 11: 462-477. 

Budyko, M. I., 1974. Climate and Life, Academic Press, 508 pp. 

Chen, T.-C., 1985. Global water vapor flux and maintenance during FGGE. 
Monthly Weather Review, 113: 1801-1819. 

Dettwiller, J. and S. A. Changnon, Jr., 1976. Possible urban effects on 

maximum daily rainfall at Paris, St. Louis and Chicago. Journal of Applied 
Meteorology, 15: 517-519. 

Eagleson, P. S., ed., 1982. Land Surface Processes in Atmospheric General 
Circulation Models. Cambridge University Press, 560 pp. 

Entekhabi, D. and P. S. Eagleson, 1989. Landsurface hydrology parameteriza- 
tion for atmospheric general circulation models: inclusion of subgrid scale 

spatial variability and screening with a simple climate model. MIT 
Department of Civil Engineering, Ralph M. Parsons Laboratory Report No. 
S25, 195 pp. 

Entekhabi, D., I. Rodriguez-Iturbe and R. L. Bras, 1991. Interannual 

variability in a stochastic model of large-scale water balance with land 
surface— atmosphere interaction, submitted to Journal of Climate. 

Fowler, W. B. and J. D. Helvey, 1974. Effect of large-scale irrigation on 
climate in the Columbia Basin. Science 184: 121-127. 

Harris, R. G., A. Thomasell, Jr. and J. G. Welsh, 1966. Studies of techniques 
for the analysis and prediction of temperature in the ocean, Part III: 
Automated analysis and prediction. Interim Report, prepared by Travelers 
Research Center, Inc. for U. S. Naval Oceanographic Office, Contract 
#N62306-1675, 97 pp. (as cited in Oort, 1983). 

Hastenrath, S. L., 1966. The flux of atmospheric water vapor over the 
Caribbean Sea and the Gulf of Mexico. Journal of Applied Meteorology, 
5: 778-788. 


130 


Henderson-Sellers, A. and K. McGuffie, 1987. A Climate Modelling Primer. 
Chichester: John Wiley &: Sons, 217 pp. 

Joussame, S., R. Sadourny and C. Vignal, 1986. Origin of precipitating water 
in a numerical simulation of the July climate. Ocean-Air Interactions , 1: 
43-56. 

Koster, R. D., J. Jouzel, R. Suozzo, G. Russell, W. Broecker, D. Read and 
P. S. Eagleson, 1986. Global sources of local precipitation as determined by 
the NASA/GISS GCM, Geophysical Research Letters, 13: 121-124. 

Koster, R. D., P. S. Eagleson and W. S. Broecker, 1988. Tracer water 
transport and subgrid precipitation variation within atmospheric general 
circulation models. MIT Department of Civil Engineering, Ralph M. Parsons 
Laboratory Report No. SI 7, 364 pp. 

Koster, R. D. and P. S. Eagleson, 1990. A one-dimensional interactive 

soil-atmosphere model for testing formulations of surface hydrology. Journal 
of Climate, 3: 593-606. 

Lettau, H., 1969. Evapotranspiration Climatonomy: 1. A new approach to 

numerical prediction of monthly evapotranspiration, runoff, and soil moisture 
storage. Monthly Weather Review, 97: 691-699. 

Lettau, H., K. Lettau, and L. C. B. Molion, 1979. Amazonia’s hydrologic cycle 
and the role of atmospheric recycling in assessing deforestation effects. 
Monthly Weather Review, 107: 227-238. 

Lufkin, D., 1959. Atmospheric water divergence and the water balance at the 
Earth’s surface. General Circulation Project Sci. Rept. No. J. 44 dd. 
Cambridge, MA: MIT Press. 

Marques, J., J. M. Santos, N. A. Villa Nova, and E. Saiati, 1977. Precipitable 
water and water vapor flux between Belem and Manaus. Acta Amazonica 7: 
355-362 (as cited in Saiati et al., 1979). 

McDonald, J., 1962. The evaporation-precipitation fallacy. Weather , 17- 
168-177. 

Oort, A. H., 1983. Global atmospheric circulation statistics, 1958-1973. NOAA 
Prof. Paper No. U, 180 pp. (+ 47 microfiches). Washington, D.C.: U. S. 

Government Printing Office. 

Paluch, I. R., 1979. The entrainment mechanism in Colorado cumuli. Journal 
of the Atmospheric Sciences, 36: 2462-2478. 

Peixoto, J. P., 1973. Atmosphere vapor flux computations for hydrological 
purposes. Report No. 20, Reports on WMO/IHD Projects, World 

Meteorological Organization (WMO), Geneva, Switzerland, 83 pp. 

Peixoto, J. P. and A. H. Oort, 1983. The atmospheric branch of the 
hydrological cycle and climate. In Variations in the Global Water Budget, 
ed. A. Street-Perrot et al., pp. 5-65. Dordrecht: D. Reidel Publ. Co. 


131 



Rasmusson, E. M., 1967. Atmospheric water vapor transport and the water 
balance of North America. Part 1. Characteristics of the water vapor flux 
Field. Monthly Weather Review, 95: 403-426. 


Rasmusson, E. M., 1968. Atmospheric water vapor transport and the water 
balance of North America. Part 2. Large scale water balance investigations. 
Monthly Weather Review, 96: 720—734. 


Rasmusson, E. M., 1971. 
using atmospheric vapor 


A study of the hydrology of eastern North America 
flux data. Monthly Weather Review, 99: 119-135. 


Rasmusson, E. M., 1977. Hydrological application of atmospheric vapor-flux 
analyses. Operational Hydrology Report No. 11, WMO, Geneva, Switzerland, 
50 pp. 


Rodriguez-Iturbe, I., D. Entekhabi and R. L. Bras, 1990. Non-linear dynamics 
of soil moisture at climate scales: stochastic analysis. Massachusetts 

Institute of Technology Center for Global Change Science Report No. 5. 


Rodriguez-Iturbe, I., D. Entekhabi and R. L. Bras, 1990. Non-linear dynamics 
of soil moisture at climate scales: chaotic analysis. Massachusetts Institute 

of Technology Center for Global Change Science Report No. 6. 

Rodriguez-Iturbe, I., D. Entekhabi and R. L. Bras, 1991. Non-linear dynamics 
of soil moisture at climate scales: 1. stochastic analysis. Water Resources 

Research, in press. 


Rosen, R. D., D. A. Salstein and J. P. Peixoto, 1979. Variability in the annual 
fields of large-scale atmospheric water vapor transport. Monthly Weather 
Review, 107: 26-27 (as cited in Oort, 1983). 

Salati E., A. Dall’Olio, E. Matsui, and J. Gat, 1979. Recycling of water in the 
Amazon Basin: an isotopic study. Water Resources Research, 15: 

1250-1258. 


Savijarvi, H. I., 1988. Global energy and moisture budgets from rawinsonde 
data. Monthly Weather Review, 116: 417—430. 

Schickedanz P. T., 1976. The effect of irrigation on precipitation in the Great 
Pl^ns. National Science Foundation Report No. NSF GI-4S871, 105 pp. 


Shiklomanov, I. A., 1989. Climate and water resources. Hydrological Sciences 
Journal, 34: 495-529. 

Shukla, J. and Y. Mintz, 1982. Influence of land-surface evapotranspiration on 
the* earth’s climate. Science, 215: 1498-1501. 

Starr, V. P. and J. P. Peixoto, 1958. On the global balance of water vapor and 
the hydrology of deserts. Tellus, 10: 189-194. 

Stidd C. K., 1967. Local moisture and precipitation. Preprint No. 45A (rev. 
July, 1968), Desert Research Institute, University of Nevada, Reno, 34 pp. 


132 


Stidd, C. K., 1975. Irrigation increases rainfall? Science 188: 279-281. 

Sutcliffe, R. C., 1956. Water balance and the general circulation of the 
atmosphere. Quarterly Journal of the Royal Meteorological Society, 82: 
385 395 . 

Taylor, G. R. and M. B. Baker, 1991. Entrainment and detrainment in 
cumulus clouds. Journal of the Atmospheric Sciences, 48: 112-121. 

Thornthwaite, C. W., 1946. The moisture-factor in climate. Transactions, 
American Geophysical Union, 27: 41 — 48. 


133 



APPENDIX A 

Computational Methods 


A.l Area integrals on the spherical surface 

We wish to approximate the integral 



(A.l) 


given values of h at nodes on a latitude longitude grid having n*m nodes 
(Figure A.l). We approximate (A.l) by a discrete sum, 

n m 

[ h(xjr) dA s l l h.. A A., w.. (A.2) 

i=l j=l 

in which h- • is the value of the function h at node (i,j), AAj. is the area of 
the rectangle with node (i,j) at its center, and w^ is a weighting term, equal 
to 1.0, 0.5, or 0.25, depending upon what fraction of the rectangle is included 
in the region. As shown in Figure A.l, 


w - w = w =w = 0.25 (corner nodes) 
w ll - w nl lm nm v 


w- = w- = 0.5, i = 1 to n 
U lm ’ 


w, . = w . = 0.5, j = 1 to m 
lj nj 


(edge nodes) (A. 3) 


w.. = 1.0 
ij 


(all other nodes) . 


The area of a rectangle on the spherical surface, AA^, is approximated as 
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AAj- 2 (A<£)(AA) cos ^ R 2 (A- 4 ) 

^ J 

in which A <p and AA are, respectively, the latitude and longitude spacing of 
the grid, in radians; ^ is the latitude at node (i,j), and R is the radius of the 
Earth, 6.37 * 10« m. 


A.2 Effective length 

The average moisture flux vector over a region, (3) has scalar 


components (Q^) and (Q^), defined as 

// Qa dA 

<V = J7d 7T 


(A.5) 


// Q\ dA 

< q a> = 71 'dr- 

where the area integrals are evaluated as described in Section A.l. A unit 


direction vector is defined 

<a> 

“ _ ll(Q)ll ' 

The angle between the positive x axis and 3 is 0, i.e. 


(A.6) 


0 = tan' 


<y 

<q a > 


(A.7) 


For purposes of this approximation, the entire region is treated as a plane 
rectangle with length Ax, where 


Ax = 


Sy 


(A-8) 


in which A is the area and Ay the north-south dimension of the region. 
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We wish to find the average length of the region in the direction of u. 
As shown in Figure A. 2, the region can be represented as n strips of width 
<5w. Each strip, i, has a length l, and the average length is 


/ _ i=n 

eff ~ 


I f i^ 


. 1 , Sv 


i =1 


(A 9) 


The unit direction vector defines at most two sides of the region as influx 
boundaries. Without loss of generality, assume that 0 < 9 < as in Figure 
A. 2. Along the north-south influx boundary, <5w = 6y cos 9 and along the 
west-east influx boundary, <5w = 6x sin 9. In the limit, taking 6w 
(consequently <5x and <5y) infinitesimal, 


^ fa — | dw = dy cos 9 + J dx sin 9 

i=l influx 0 n 


influx 0 
boundary 


(A. 10) 


Noting that 16* = A, (A. 9) becomes 


I _ A 

eff Ay cos0 + Ax sin 0 ‘ 


(A. 11) 


More generally, to allow for u in other quadrants, 




Ay | cos 9 1 + Ax | sin 9 \ 


(A. 12) 
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A.3 Effective moisture influx 

To determine (wu)^ we resolve <3 along the influx boundary into a 
component parallel to u and a component perpendicular to u. The parallel 
component is <3-u. Following the same argument as above, the length 
weighted average <3-u is 

n 

y (3-u). L6 w 

= H ' ' (A.13) 


In the limit, integrating along the influx boundary, 



(A.14) 


in which 7 is the distance along the influx boundary, 7^ and 7^ are the 
endpoints of the influx boundary, (I7 is the unit vector defining the direction 
in which the path is traced, and a is the angle between (I7 and u. As shown 


in Figure A.3, a = 6 along the west-east boundary, and a = Q + ^ along the 
north-south boundary, when 0 < 0 < ^. 
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APPENDIX B 


Computer Programs 
Programs include: 


1. AREA. Computes the area of a region, given latitude, longitude, and 

weight factor for each node. 

2. BFLUXW. Computes the net atmospheric vapor flux over a specified 

boundary on the plane 

3. QEFF. Computes areally-averaged moisture flux vector components. 

4. WU. Computes (wu) eff , Equation (2.16). 


Common to all the computer routines is a database of the components of <3 in 
73 x 73 (ILAT.JLON) arrays, where ILAT is the latitude index (1 at 90° S to 
73 at 90° N) and JLON is the longitude index (1 at 180° E to 73 at 180° W). 
The components were computed by trapezoidal rule integration of Equation 


141 



area . f 


Wed May 22 13:25:32 1991 


1 


PROGRAM AREA 
C 

C Computes area, given latitude, longitude, 

C and included part (SQUSE) of grid square surrounding each node. 

C (because some nodes are at corners & do not contribute whole elements) 
C Reads an input file: N, M, beginning indices (GFDL Grid 
C start top left) and N row by M column matrix of area fractions. 

C units: area fractions are in eigths. (0, 1, 2, 3, 4 . . . 8) 

C Computes an estimate of the total area defined by the grid. 

C Note: "estimate" means that the area (dy)x(dx) is estimated as 
C dA = (R*dPhi) (R*dLam*cos (PhiNode) ) , PhiNode is the latitude at 
C the node, i.e. the center of the grid element. 

C 

C Area is in m~2. 

C 

CHARACTER* 60 TITLE 
CHARACTER* 5 FNSET 
CHARACTER* 18 FNSQAR, FNAREA 
REAL PI , PHI , SQUSE (73, 73) 

PARAMETER (AR=6.371E6, PI=3. 14159) 

DATA DPHI , DLAM / 2.5, -5./ 

DATA LAT2 / 37 / 

C 

C Convert to radians 

DPHI=DPHI / 1 8 0 . *PI 
DLAM=DLAM/ 180. *PI 
C 

WRITE (6,*) 'Name of dataset?' 

READ (5,*) FNSET 
FNSQAR=FNSET/ / ' . sq' 

FNAREA=FNSET//' . ar' 

OPEN (UNIT=12, FILE=FNSQAR, STATUS"' OLD' ) 

READ (12, 100) TITLE 
READ (12,*) N, M, LATB, LONE 
LATF=LATB-N+ 1 
LCNF=LONB-M+l 
SUM=0 . 0 

DO 10 ILAT=LATB, LATF, -1 

PHI=REAL ( I LAT-LATZ ) *DPHI 
COSPHI=COS (PHI) 

WRITE (6,*) PHI , COSPHI 

READ (12, * ) (SQUSE (ILAT, JLON) , JLON=LONB, LONF , -1) 

WRITE (6,*) (SQUSE (ILAT, JLON) , JLON=LONB, LONF, -1 ) 

DO 20 JLON=LONB, LONF, -1 

TERM=COSPHI* SQUSE (ILAT, JLON) /8 . 

SUM=SUM+TERM 
20 CONTINUE 

10 CONTINUE 

CLOSE (12, STATUS"' KEEP' ) 

C 

RAREA= SUM* AR*AR* DPHI *ABS (DLAM) 

OPEN (UNIT=13, FILE=FNAREA, STATUS"' NEW' ) 

WRITE (13, 500) FNSQAR 
WRITE (13, 501) RAREA 
STOP 
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100 FORMAT { A 6 0 ) 

500 FORMAT (IX, 'Area Computation for ' , A1 8 , 3X, ' (m~2 ) ' / ) 

501 FORMAT (IX, Ell. 4) 

END 
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PROGRAM BFLUXW 
Q 

c Computes the net atmospheric vapor flux over a specified boundary 
C on the plane (spherical), using GFDL' s 73x73 (Lat,Lon) grid. 

C Created 7/30/90, KB. Modified for unix fortran 2/28/91. 

C Applies bivariate interpolation (allows diagonal boundary segments) . 

C The two components of the vector must be specified at each node 
C along the boundary, and at all corner nodes of a grid when the 
C boundary segment is diagonal. 

C The net INflux over an infinitesimal segment of the boundary is 
C equal to the dot product of the flux vector and the inward normal , 

C vector; in this case, with linear boundary segments, each segment s 
C contribution is the dot product of inward normal and average flux 
C vector for the segment . 

C The resulting total over the whole boundary is -divV x A, 

C where V is the vector quantity and A is the area of the region. 

C Influx and Cutflux are both computed, 

c 

C For the variables, the final index distinguishes between the 
C value of the quantity, and its variance (sigma~2), as computed by 
C standard propagation of error. The quantity is then expressed 
C as +- 2* sigma. 

C 

C Components of the vector moisture flux (QPHI, QLAM) have units 
C g m" ( -1 ) s" ( - 1 ) 

C Distances (DELX , DELY , BNX , ENY , BX , BY ) have units m 
C 

CHARACTER* 6 0 TITLE, CJUNK 
CHARACTER* 5 FNSET 

CHARACTER* 1 8 FNBNOD, FNAREA, FNSIGF, FNTXT , FNSGFX 
CHARACTER* 3 FNAME (17) 

CHARACTER* 30 FNLAM, FNPHI 
INTEGER LAT(IOO), LCN(lOO) 

REAL QPHI (73, 73, 2) , QLAM ( 7 3 , 7 3 , 2 ) , FX(2), FY(2), QIN(2), 

& QNET (2 ) , SUMIN (2 ) , SUMOUT(2), 

& BNX (100), ENY (100), SFIN (17), SFOUT(17), SEGFLX (100, i 7 ) , 

& SQNET (17) 

REAL PHI 1 , PHI 2 

PARAMETER (AR=6.371E6, PI=3. 14159, RHOW=1000 . , SECDAY=864 00 . ) 

DATA DPHI , DLAM / 2.5, -5./ 

DATA FNAME / ' jan' , ' feb' , ' mar ' , ' apr ' , ' may' , ' jun' , 

& ' j.ul ' , ' aug' , ' sep' , ' oct' , ' nov' , ' dec' , 

& 'djf' , 'mam' , ' j ja' , ' son' , 'yer' / 

C 

100 FORMAT (A60 ) 

101 FORMAT (' 1' ,A60, //, IX, 'Month' , 13, IX, '(', A3, ')' //, , ^ 

&1 6X, ' Flux, kg rrr-1 s~-l Scaled inward Segment , flux' , / 

&16X, ' normal vector, m kg s~-l' ,// 

&' Node K K+l F (x) F (y ) B(x) B(y) Qin 

& IX 7 2 ( ' = 9 ) ) 

102 FORMAT ( IX , 14 , 1 6 , 4X, 5 (El 0 . 3 , IX) , ' +-' , IX, El 0 . 3 ) 

105 FORMAT (7E10.3 ) 

106 FORMAT ( / , IX, ' Area = ',E10.3,' m~2',// 

&1X, ' <Qin> = ' , El 0 . 3 , ' mm/day IX, E9 . 3, 

& 12X, ' Qin = ' , E 10 . 3, IX, ' +-' , IX, E10 . 3, /, 
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&lX,'<Qout> = ',E10.3,' mm/day IX, E9 . 3 , 

& 12X, ' Qout = ' ,E10. 3, IX, , 1X,E10.3, //, 

&lX,'<Qnet> = ' , E10 . 3, ' mm/day IX, E9 . 3, 

& 12X,'Qnet = ' , E10 . 3, IX, ' +-' , IX, E10 . 3, /) 

600 FORMAT (IX, ' Boundary Flux Computation for ' , A18, 3X, ' (kg s~-l)'/') 

601 FORMAT (IX, 12E10 . 3) 

606 FORMAT (IX, 4 (E10 . 3,20X) ) 

602 FORMAT (/, IX, ' Boundary Flux Computation for ',A18,3X, 

& ' (mm mo "'-l ) ' / ) 

603 FORMAT ( IX, 12F1 0 . 1 ) 

607 FORMAT ( IX , 4 ( FI 0 . 1 , 2 OX ) ) 

604 FORMAT (17 (IX, F5. 1) ) 

605 FORMAT (/) 

500 FORMAT (A60,/) 

501 FORMAT (IX, Ell. 4) 

C 

C Additional conversion factor 
SECMON=3 0 . 4 * SECDAY 
C Convert to radians 

DPHI=DPHI/180 . "PI 
DLAM=DLAM/1 80 . *?I 

Input job data. FNBNOD has boundary node specifications. 

FNAREA has region's area. 

WRITE (6,*) 'Name of Data Set?' 

READ (5,*) FNSET 
FNBNCD=FNSET//' .bn' 

FNAREA=FNSET/'/' .ax' 

FNTXT=FNSET//' .txt' 

FNSIGF=FNSET/ / ' .fix' 

FNSGFX=FNSET//' . seg' 

OPEN (UNIT=8 , F I LE= FNBNOD, STATUS=' OLD' ) 

READ (8,100) TITLE 
N=0 

10 IF (N .GE. 999) GOTO 1 1 
NPTS=N 

READ (8,*) N, LAT(N), LON (N) 

GOTO 10 

11 CONTINUE 
CLOSE (8, STATUS 55 ' KEEP ' ) 

OPEN (UNIT=8 , FILE=FNAREA, STATUS='OLD' ) 

READ (8,500) C JUNK 
READ (8,501) SAREA 
CLOSE (8, STATUS=' KEEP' ) 

== BOUNDARY LOOP ================= ==== = = ============================ , 

C 

DO 60 K=1 , NPTS-1 

C Get delta x, delta y for bdy segment K to K+l . 

KP 1 =K+ 1 
I1=LAT (K) 

I2=LAT ( KP 1 ) 

Jl=LON (K) 

J2=LON ( KP 1 ) 

PHI 1=DPHI * (11-37) 
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PHI2=DPHI * (12-37) 

PHI AV= (PHI 1 + PHI2 ) / 2 . 

DELX= ( J2-J1 ) *DLAM*AR*COS (PHIAV) 

DELY= (12-11) *DPHI * AR 

C Components of the inward normal vector. 

BNX (K) =-DELY 
BNY (K) = DELX 
60 CONTINUE 
C 

C == End boundary loop ================================== == 

C 

OPEN (UNIT=6 , FILE=FNTXT, STATUS=' NEW' ) 

== MONTH LOOP ======================================== ==:== == ==:=====!:= 

DO 70 IMON-1,12 
C Load component arrays. 

FNLAM=FNAME (IMON) / / ' _qlam. dat ' 

FNPHI=FNAME (IMON) II' _qphi . dat' 

OPEN (UNIT=9 , FILE=FNLAM, STATUS*' OLD' ) 

OPEN (UNIT=10, FILE=FNPHI , STATUS*' OLD' ) 

DO 75, K-1,2 
DO 751=1, 73 

READ (9,105) (QLAM ( I , J, K) , J=l,73) 

READ (10, 105) (QPHI (I, J,K) ,0=1,73) 

75 CONTINUE 

CLOSE (9, STATUS*' KEEP' ) 

CLOSE (10, STATUS*' KEEP' ) 

C 

SUMIN (1 ) =0 . 

SUMIN (2) =0 . 

SUMOUT (1 ) =0 . 

SUMOUT (2) =0 . 

WRITE (8,101) TITLE, IMON, -NAME (IMON) 

C Compute and record influx for each bdy segment. 

DO 71 K=1 , NPTS-1 
KP1=K+1 
I 1=LAT (K) 

I2=LAT (KP1 ) 

Jl=LON (K) 

J2=LON (KP1 ) 

C Averaged flux vector components, divide by 1000 to convert g to Kg 
FX ( 1 ) * (QLAM (II, J1 , 1 ) +QLAM (12, J2, 1) ) / 3. + 

& (QLAM (II , J2, 1) +QLAM (12, J1 , 1 ) ) / 6. 

FX (2 ) * (QLAM (II , J1 , 2) +QLAM (12, J2, 2) ) / 9. + 

& (QLAM ( I 1 , J2 , 2 ) +QLAM ( 12 , J1 , 2 ) ) / 36. 

FX ( 1 ) * FX ( 1 ) /1000 . 

FX (2 ) = FX (2) /1000 . /1000 . 

FY ( 1 ) = (QPHI (II, Jl, 1) +QPKI (12, J2, 1) ) / 3. + 

& (QPHI (II, J2, 1) +QPHI (12, Jl, 1) ) / 6. 

FY (2) = (QPHI (II, J1,2)+QPHI (12, J2,2) ) / 9. + 

& (QPHI (II, J2, 2) +QPHI (12, Jl, 2) ) / 36. 

FY (1 ) = FY (1) /1C00 . 

FY (2 ) = FY (2) /1000./100C. 

C Dot product . Sums . 
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BX=BNX (K) 

BY=BNY (K) 

QIN ( 1 ) =BX *FX ( 1 ) + BY* FY ( 1 ) 

QIN (2) =BX*BX*FX (2) +BY*BY*FY (2) 

IF (QIN (1 ) . GE . 0 . ) THEN 

SUM IN ( 1 ) =SUMIN ( 1 ) +QIN ( 1 ) 

SUMIN (2) =SUMIN (2) +QIN (2) 

ELSE 

SUMOUT ( 1 ) =SUMOUT ( 1 ) +QIN ( 1 ) 

SUMOUT (2) =SUMOUT (2) +QIN (2) 

END IF 

C Output 

WRITE (8,102) K, KP1, FX (1) , FY (1) , BX, BY, QIN (1) , 2 . *SQRT (QIN (2) ) 
SEGFLX (K, IMON) =QIN (1) /1E+6 
71 CONTINUE 

Unit conversion . m to mm,* kg water to in^3, sec to days. 

Divide by area to get <Q> . 

UCONV=l 000 . /RHOW*SECDAY 
RFACT=UCONV / SAREA 
QNET ( 1 ) =SUMIN (1 ) + SUMOUT (1 ) 

QNET (2) =SUMIN (2) + SUMOUT (2) 

SUMOUT ( 1 ) =ABS ( SUMOUT ( 1 ) ) 

WRITE (8,106) SAREA, 

& SUMIN ( 1 ) *RFACT, 2 . * 5QRT ( SUMIN (2 ) ) *RFACT, 

& SUMIN (1), 2 . *SQRT (SUMIN (2) ) , 

& SUMOUT (1) *RFACT, 2 . * SQRT ( SUMOUT ( 2 ) ) *RFACT, 

& SUMOUT (1 ), 2 . *SQRT (SUMOUT (2) ) , 

& QNET (1) *RFACT, 2 . * SQRT (QNET (2 ) ) *RFACT, 

& QNET ( 1 ) , 2. * SQRT (QNET (2) ) 

SFIN (IMON) =5UMIN ( 1 ) 

SFOUT (IMON) = SUMOUT (1) 

SQNET ( IMON) =QNET ( 1 ) 

70 CONTINUE 

C = = END MONTH LOOP ====»-=========================== ================== 

C 

CLOSE (8, STATUS^' KEEP' ) 

r* 

UCQNV=1000 . / RHOW* SECMON 
RFACT=UCONV/ SAREA 

c Write to data file for balance ccir.ps . 

OPEN (UNIT=10 , FILE=FNSIGF, STATUS=' NEW' ) 

WRITE (10,600) FNBNOD 

WRITE (10,601) (SFIN(I) , 1=1, 12) 

WRITE (10,601) (SFOUT (I) , 1=1, 12) 

WRITE (10,601) (SQNET (I) , 1 = 1, 12) 

WRITE (10,605) 

WRITE (10,606) (SFIN (I ), 1=13, 16) 

WRITE (10,606) (SFOUT(I) , 1=13, 16) 

WRITE (10,606) (SQNET (I) , 1=13, 16) 

WRITE (10,605) 

WRITE (10,606) SFIN (17) 

WRITE (10,606) SFOUT (17) 

WRITE (10,606) SQNET (17) 

WRITE (10,602) FNBNOD 
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WRITE (10,603) (SFIN (I) *RFACT, 1=1, 12) 

WRITE (10,603) ( SFOUT ( I ) *RFACT , 1=1 , 12 ) 

WRITE (10,603) (SQNET (I) *RFACT, 1=1, 12) 

WRITE (10,605) 

WRITE (10,607) (SFIN (I) *RFACT, 1=13, 16) 

WRITE (10,607) (SFOUT (I) *RFACT, 1=13, 16) 

WRITE (10,607) (SQNET (I) *RFACT, 1=13, 16) 

WRITE (10,605) 

WRITE (10,607) SFIN (17) ‘REACT 
WRITE (10,607) SFOUT (17) *RFACT 
WRITE (10,607) SQNET (17) *RF ACT 
CLOSE (10, STATUS=' KEEP' ) . 

C Write to segment flux file for barplotting. 

OPEN (UNIT=10, FILE=FNSGFX, STATUS=' new' ) 

DO 65 K=1 , NPTS-1 

WRITE (10,604) (SEGFLX (K, I MON) , IMON=l, 17) 
65 CONTINUE 

CLOSE (10, STATUS-' keep' ) 

C 

STOP 

C 

END 
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PROGRAM QEFF 
C 

C 3/12/91 KB 

C Given latitude, longitude at nodes; 

C and included part (SQUSE) of grid square surrounding each node. 

C {because some nodes are at corners & do not contribute whole elements) , 

C computes areally-averaged values of the moisture flux vector components. 
C Output: <Qlam>, <Qphi>, | <Q> | , and theta. 

C Reads an input file: N, M, beginning indices (GFDL Grid — 

C start top left) and N row by M column matrix of area fractions. 

C units: area fractions are in eigths. < 0 , 1 , 2 , 3 , 4 . . . 8 ) 

C Computes an estimate of the total area defined by the grid. 

C Note: "estimate" means that the area (dy)x(dx) is estimated as 
C dA = (R*dPhi) (R*dLam*cos (PhiNode) ) , PhiNode is the latitude at 
C the node, i.e. the center of the grid element. 

C 

C Area is in m^2. 


C Addition 3/29/91: 

7 Computes average components and theta for each GFDL 

C in the region. (For checking the validity of assumi 


c 


c 


CHARACTER* 60 TITLE 
CHARACTER* 5 FNSET 


CHARACTER* 18 FNSQAR, FNCUT, FNEOUT, FNLAM, FNPHI 
CHARACTER*3 FN ( 12 ) 

REAL PI , PHI , SQUSE (73,73) , TERM (73,73), ETERM(73,73) 
& QLAM (73, 73, 2) , QPHI (73, 73, 2) 

PARAMETER (PI=3. 14159) 

DATA DP HI , DLAM / 2.5, -5./ 

DATA LATZ / 37 / 

DATA FN /' jan' , ' feb' , 'mar' , ' apr' , ' may' , ' jun' , 


' jul ' , ' aug' , ' sep' , ' oct' , ' nov' , ' ae 


Convert to radians 

DPHI=DPHI/180 . *PI 
DLAM=DLAM/ 180. *PI 


WRITE (6,*) ' >> Name of dataset?' 

READ (5,*) FNSET 
FNSQAR=FNSET/ / ' . sq' 

FNOUT= FNSET//' .qef' 

FNEOUT= FNSET//' .eleqef' 

OPEN (UNIT=12 , FILE=FNSQAR, STATUS*' OLD' ) 

READ (12, 100) TITLE 

READ (12,*) N, M, LATB, LONB 

LATF=LATB-N+1 

LONF=LONB-M+l 

WASUM=0 . 0 

DO 10 ILAT=LATB, LATF, -1 

PHI=REAL (ILAT-LATZ) *DPHI 
COSPHI=COS (PHI) 

READ (12, * ) (SQUSE (ILAT, JLON) , JLON=LONB, LONF, -1) 
DO 20 JLON=LONB, LONF, -1 

TERM (ILAT, JLON) =COSPHI * SQUSE ( ILAT, JLON) /8 . 


lement 
f c rmi ty ) 
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ETERM (ILAT, JLON) =COSPHI 
WASUM=WASUM+TERM ( ILAT , JLCN) 

20 CONTINUE 

ONTINUE 

LOSE (12, STATUS*' KEEP' ) 

OPEN (UNIT=13, FILE=FNO T JT, STATUS = 'NEW' ) 

OPEN (UNIT=15, FILE=FNEOUT, STATUS*' NEW' ) 

WRITE (13,200) FNSET 
DO 50 IMON=l , 12 

FNLAM=FN (IMON) / / ' _qlam . dat ' 

FNPHI=FN (IMON) / / ' _qphi . dat ' 

OPEN (Unit=l 4 , f ile=FNLAM, status*' old' ) 

OPEN (Unit=16, f ile=FNPHI , status*' old' ) 

DO 60 K=1 , 2 
DO 60 1=1,73 

READ (14,105) (QLAM(I, J,K) , J-1,73) 

READ (16,105) (QPHI (I, J,K) , J=l,73) 

60 CONTINUE 

CLOSE (14) 

CLOSE (16) 

QLAMSM=0 . 0 
QPHI SM=0 . 0 

DO 65 ILAT=LATB, LATF, -1 

DO 65 JLON=LONB, LONF, -1 

QLAMSM=QLAMSM+ TERM (ILAT, JLCN) *QLAM (ILAT, JLON, 1 ) 
QPHISM=QFHISM+TERM (ILAT, JLON) *QPHI (ILAT, JLON, 1 ) 
65 CONTINUE 

QLAMAV=QLAMSM/WASUM 
QPHIAV=QPHI SM/WASUM 

QAV=SQRT ( QLAMAV *QLAMAV+QPKIAV* QPHIAV) 

THETA=ATAN (QPHI AV 1 QLAMAV ) *180 . /PI 
C ATAN returns theta in range (-90,90) 

IF (THETA. LT. 0 .AND. QLAMAV. LT. 0) THETA=THETA+1 8 0 . 

IF (THETA. GT. 0 .AND. QPHIAV. LT. 0) THETA=THETA-1 8 0 . 
WRITE (13,500 IMON, QLAMAV, QPHIAV, QAV , THETA 

Q 

C eleir.ent-wise averaging 

C 

WRITE (15, 700) IMON 
DO 75 ILAT=LATB, LATF+1, -1 
DO 75 JLON=LONB, LONF+1 , -1 
I1=ILAT 
I2=ILAT-1 
Jl= JLON 
J2=JLON-l 


& 

EWSUM=ETERM ( 11 , J1 ) +ETERM ( 11 , J2 ) + 
ETERM (12, Jl) +ETERM(I2, J2) 



ESUMPHI=ETERM ( 11 , Jl ) 

*QPHI (11 , Jl, 1 

) + 

& 

ETERM ( I 1 , J2 ) 

*QPHI (11, J2, 1 

) + 

& 

ETERM(I2, Jl) 

*QPHI (12, Jl, 1 

) + 

& 

ETERM (12, J2) 

*QPHI (12, J2, 1 

) 


ESUMLAM=ETERM ( 1 1 , Jl ) 

*QLAM(I1, Jl, 1 

) + 

& 

ETERM (11, J2) 

*QLAM (11, J2, 1 

) + 

& 

ETERM(I2, Jl) 

*QLAM(I2, Jl, 1 

) + 
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& ETERM(I2, J2) *QLAM(I2, J2, 1) 

EAVPHI=ESUMPHI /EWSUM 
EAVLAM=ESUMLAM/EWSUM 

WRITE (15,705) ILAT, JLON, EAVLAM, EAVPHI 
75 CONTINUE 

C 

50 CONTINUE 

CLOSE (13) 

CLOSE (15) 

STOP 

C 

100 FORMAT (A60) 

105 FORMAT (7E10.3) 

200 FORMAT (IX, 'Area Mean Q Vector for ',A18,//, 

& IX, 'Month <Q (lam) > <Q(phi)> I <Q> I 

& ' t g / s/m] 

500 FORMAT (IX, I5,2X, 3(2X,E10.3) ,F9.1) 

501 FORMAT (IX, Ell. 4) 

'00 FORMAT ( / , IX , ' IMON=' , 1 3 , 2X, ' Element Averages',//, 
& IX, 'Top LH Node <Q ( 1 am) > <Q(phi)>') 

^705 FORMAT ( IX , 1 2 , IX, 12 , 5X , 2 ( 2X , El 0 . 3 ) ) 

END 


Theta' , / 
[deg] ' , // ) 
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PROGRAM WU 
C 

c Computes the "wu" term for Budyko recyling analysis, 
c 

C Here, wu = a weighted average of the inward moisture flux vector 
C component in the direction of the average moisture flux vector over 
C the region. 

C 

C The unit direction vector is U, with components ULAM, UPHI 
C 

C Qphi , Qlambda , and Qin, parallel (wu) have units g mM-l) s"(-D 
c 

C This routine works only for rectangular regions. 

Q 

CHARACTER* 18 FNAR, FNQEF, FNSQ, FNVJU, FNBN, FNCHK, FNLAM, FNPHI 
CHARACTER* 5 FNSET 
CHARACTER* 3 FNAME ( 12 ) 

CHARACTER* 60 TITLE 
INTEGER IK (73) , JK (73) 

REAL QLAM (73,73) , QPHI ( 73 , 7 3 ) , RL ( 7 3 ) , QP ( 7 3 ) , DW ( 7 3 ) 

REAL DELY, YREL, YMAX, XMAX, DELX 
PARAMETER ( RPI=3 . 1 4 1 59 , AR- 6.37E+6 ) 

DATA FNAME / ' jan' , ' f eb' , ' mar' , ' apr ' , ' may ' , ' jun' , 

4 ' jul ' , ' aug' , ' sep' , ' cct' , ' nov' , ' dec' / 


1C5 FORMAT (7E10.3) 

200 FORMAT (//, IX, Ell . 4 ) , 

300 FORMAT ( ' RAREA= ' , El 0 . 3 , 2X , ' YMAX= ' , El 0 . 3 , 2X, ' XMAX -',E10.o,/ 
& 1 9X , ' DELY= ' ,E10.3,2X, 'DELX =',E10.3) 

305 FORMAT (/,'** IMON= ',14,' **', 

4 //,'QLAMAV= ' , E10 . 3, 2X, ' QPHIAV= ',E10.3,2X, 

4 ' QAV= ' ,E10.3,2X, 'THETA= ',F6.1) 

310 FORMAT (/,'ABCOS= ' , F8 . 4 , 2X, ' ABSIN= ',F8.4) 

500 FORMAT (/////) 

505 FORMAT ( 8X , 3 (2X , El 0 . 3) , F9 . 1 ) 

600 FORMAT ( ' ULAM=' ,F6.2,3X,'UPHI=' ,F6.2) 

601 FORMAT ( / , ' ILATB= ' , 1 3 , 2X , ' I INC=' , 1 3 , 2X, ' ILATF=' , 1 3 , 4X , 

4 ' JLONB=' , 13, 2X, ' JINC=' , 13, 2X, ' JLONF=' , 13) 

602 FORMAT ( / , ' X-segment integration',/, 

4 ' K JK DELW RLZ RL1 AREA 

4 WTZ WT1 QPBAR' ) 

603 FORMAT (12 , 14 , 6 (E10 . 3) , 2 (F6 . 3) , Ell . 3) 

604 FORMAT (/,'CRIT VAL = ' , E 9 . 3 , 2X , ' RLMAX= ',E9.3) 

605 FORMAT ( / , ' Y-segment integration',/, 

4 ' K IK DW (K) RLZ RL1 AREA 

4 WTZ WT1 QPBAR' ) 

607 FORMAT ( / , 31X , ' ASUM= ' , E9 . 3 , 27X, ' QPSUM= ',E9.3,/ 

&/, 31X, ' IMON= ' , 14 , 6X, ' ** QPAVG= ',E9.3,' **') 

700 FORMAT ( IX, ' Weighted Mean Q (Parallel )' , 

4 3X/[g m A (-l) s~(-l)]',3X,A18,//, 

4 IX , ' IMON <Qp>' ) 

705 FORMAT ( IX, 13, 2X, E9 . 3) 

C 

WRITE (6,*) ' >> Name of dataset?' 

READ (5,*) FNSET 


QPZ 


QPZ 


QP1 


QP1 
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FNAR = FNSET//' . ar' 

FNQEF= FNSET//' .qef' 

FNSQ = FNSET//' . sq' 

FNWU = FNSET//' .wu' 

FNBN = FNSET//' .bn' 

FNCHK= FNSET//' . cnk' 

OPEN (unit=25, file=FNCHK, status=' unknown' ) 

OP EN (unit=27 , fi 1 e=FNWU, st a t us= ' new' ) 

WRITE (27,700) FNSET 

OP EN ( unit = l 0 , f ile=FNSQ, status= ' old' ) 

READ(10,*) TITLE 

READ (10,*) N, M, LATB, LONB 

LATF=LATB-N+1 

LONF=LONB-M+l 

CLOSE (10) 

OPEN (unit=10, f ile=FNAR, status=' old' ) 

READ (10, 200) RAREA 
CLOSE (10) 

C 

C Planar distances. Xmax is the average longitudinal dimension 
C of the spherical-rectangular region. Xmax* Ymax=Area . 

C 

DELY=AR*RPI / 1 8 C . *2.5 
YMAX=REAL (N-l ) *DELY 
XMAX=RAREA/YMAX 
DELX=XMAX/REAL (M-l) 

WRITE (25,300) RAREA, YMAX, XMAX, DELY, DELX 
C 

OPEN (unit=12, file=FNQEF, status=' old' ) 

READ (12,500) 

r* 

C = = MONTH LOOP ======================================= == == ====; 

C 

DO 100 IMON=l ,12 

c 

C Read the Q vector comDonents for this month. 

C 

FNLAM=FNAME (IMON) //'_qlam. dat' 

FNPHI=FNAME ( IMON) / / ' _qphi . dat ' 

OPEN (UNIT=2 0 , FI LE=FNLAM, STATUS = ' old' ) 

OPEN (UNIT=22 , FILE=FNPHI, STATUS=' old' ) 

DO 75, 1=1,73 

READ (20,105) (QLAM (I, J) , J=l, 73) 

READ (22,105) (QPHI (I, J) , J=l, 73) 

75 CONTINUE 

CLOSE (20) 

CLOSE (22) 

C 

C Read the region-specific direction data for this month. 

C Compute components of the unit direction vector. 

C 

READ (12, 505) QLAMAV, QPHIAV, QAV, THETA 
write (25,305) IMON, QLAMAV, QPHIAV, QAV, THETA 
ULAM=QLAMAV/QAV 
UPHI=QPHIAV/QAV 
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WRITE (25, 600) ULAM, UPHI 
C 

C Determine direction for integration along inward boundary. 

C ( JLON, ILAT ) is the x,y corner at which to start. 

C 

If (ULAM.GT. 0 . ) THEN 
JLONB=LONF 
JLONF=LONB 
JINC-1 
Else 

JLONB-LONB 

JLONF=LONF 

JINO-1 

Endif 

I F (UPHI . GT . 0 . ) THEN 
I LATB=LATF 
I LATF=LATB 
IINC-1 
ELSE 

I LATB=LATB 
I LATF=LATF 
I INC=- 1 
ENDIF 

WRITE (25,601) ILATB , IINC, ILATF , JLONB , JINC, JLONF 
C 

C == Integrate along the X-directicn segment first. 

C 

C Initialize sums 
C 

ASUM=0 . 0 
QASUM=0 . 0 
C 

C If THETA is 180 or -180, skip the x-direction segment. 

C 

THETA=THETA*RPI/180 . 

ABCOS=ABS (COS (THETA) ) 

ABSIN=ABS (SIN (THETA) ) 

WRITE (25,310) ABCOS , ABSIN 
IF (ABS (THETA) .eq. 180. ) GOTO 150 

Compute Xcrit and fit it into the sequence of nodes for integration. 
MCOUNT=M 

XCRIT=YMAX/ABS (TAN (THETA) ) 

IF (XCRIT . LT . XMAX) MCOUNT=M+l 
RLMAX=YMAX/ABSIN 
WRITE (25,604) XCRIT, RLMAX 
IFLAG=0 
K-l 

KINT=0 
JK (1) = JLONB 
XREL=0 . 

RL (1) =0 . 

QP (1) =QLAM (ILATB, JK (1) ) *ULAM + QPHI (ILATB, JK (1 )) *UPHI 
120 IF (K.eq.MCOUNT) GOTO 121 
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KINT=KINT+1 

C XREL=XREL+DELX 

XREL=REAL (KINT) *DELX 
IF (XREL.LT.XCRIT) THEN 
K=K+ 1 

JK (K) =JK (K-l ) +JINC 
DW (K) =DELX*ABSIN 
RL (K) =XREL/ABCOS 

QP (K) =QLAM ( ILATB, JK (K) ) *ULAM+QPHI (ILATB, JK (K) ) *UPHI 

ELSE 

IF (IFLAG.LT. 1) THEN 
K=K+1 
JK (K) =99 
RL (K) =RLMAX 
XDIF=XCRIT-XREL+DELX 
DW (K) =XDIF*ABSIN 
K=K+1 

JK (K) =JK (K-2 ) +JINC 
RL (K) =RLMAX 
XDIF=XREL-XCRIT 
DW (K) =XDIF*ABSIN 

QP (K) =QLAM (ILATB, JK (K) ) *ULAM+QPHI (ILATB, JK (K) ) *UPHI 
QP (K-l ) =QP (K-2) * (XDIF/DELX) +QP (K) * ( 1 -XDIF/DELX ) 
IFLAG=1 
ELSE 

K=K+1 

JK (K) = JK (K-1J+JINC 
RL (K) =RLMAX 
DW (K) =DELX*ABSIN 

QP (K) =QLAM (ILATB, JK (K) ) *ULAM+QPHI (ILATB, JK (K) ) *UPHI 
ENDIF 
END IF 
GOTO 120 
121 CONTINUE 

n 

WRITE (25,602) 

DO 110 K=2 , MCOUNT 
RLZ=RL (K-l) 

RL1=RL (K) 

QPZ=QP (K-l) 

QP1=QP (K) 

AREA = (RLZ+RL1 ) *DW (K) /2 . 

DENOM = 3* (RLZ+RL1 ) 

WTZ = (2*RLZ+RL1 ) /DENOM 
WT1 = (RLZ+2*RL1) /DENOM 
QPBAR = WTZ*QPZ + WT1*QP1 
ASUM=ASUM+AREA 
QASUM=QASUM+AREA* QPBAR 

WRITE ( 25 , 603) K , JK ( K ) ,DW(K) , RLZ, RL1 , AREA, QPZ , QP1 , WTZ , WT1 , 

& QPBAR 

110 CONTINUE 

150 CONTINUE 

C 

== The Y-direction segment. 
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Compute Ycrit and fit it into the sequence of nodes for integration. 
NC0UNT=N+1 

YCRIT=YMAX-XMAX*ABS (TAN (THETA) ) 

RLMAX-XMAX/ABCOS 

IFLAG=0 

IF (YCRIT. LE. 0 . ) THEN 
YCRIT=0 . 

RLMAX=YMAX/ABSIN 
NCOUNT=N 
IFLAG=1 
END IF 

WRITE (25,604) YCRIT, RLMAX 
K=1 

KINT=0 

IK ( 1 ) =ILATB 

YREL-0 . 

RL (1 ) =RLMAX 

QP (1) =QLAM (IK (1) , JLONF) *ULAM + QPHI (IK (1 ) , JLONF) *UPHI 
130 IF (K.eq.NCOUNT) GOTO 131 

KINT=KINT+1 

C YREL-YREL+DELY 

YREL=REAL (KINT) *DELY 
IF (YREL. LT. YCRIT) THEN 
K-K+l 

IK (K) -IK (K-l) +IINC 
DW (K) =DELY*ABCOS 
RL (K) =RLMAX 

QP (K) -QLAM(IK(K) , JLONF ) ’ULAM+QPHI (IK(K) , JLONF) *UPHI 

ELSE 

IF (IFLAG.LT. 1) THEN 
K=K+ 1 
IK (K) =99 
RL (K) =RLMAX 
YDIF=YCRIT-YREL+DELY 
DW (K) =YDIF*ABCOS 
K=K+1 

IK (K) —IK (K-2) +IINC 
RL (K) = ( YMAX-YREL) /ABSIN 
YDIF-YREL- YCRIT 
DW (K) =YDIF*ABCOS 

QP (K) -QLAM (IK (K) , JLONF) *ULAM+QPHI (IK (K) , JLONF) *UPHI 
QP (K-l) =QP (K-2) * (YDIF/DELY) +QP <K) * (1-YDIF/DELY) 
IFLAG-1 
ELSE 

K-K+l 

IK (K) —IK (K-l ) +IINC 
RL (K) = (YMAX-YREL) /ABSIN 
DW (K) =DELY*ABCOS 

QP (K)-QLAM(IK(K) , JLONF) *ULAM+QPHI (IK(K) , JLONF) *UPHI 
END IF 
END IF 
GOTO 130 
131 CONTINUE 

C 
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C Y-segment integration Loop. 

C 

WRITE (25,605) 

DO 140 K=2 , NCOUNT 
RLZ=RL (K-l) 

RL1=RL (K) 

QPZ=QP (K-l) 

QP1=QP (K) 

AREA = (RLZ+RL1) *DW(K) /2. 

DENOM = 3* (RLZ+RL1 ) 

WTZ = (2*RLZ+RL1) /DENOM 
WT1 = (RLZ+2*RL1 ) /DENOM 
QPBAR = WTZ*QPZ + WT1*QP1 
ASUM=ASUM+AREA 
QASUM=QASUM+ AREA* QPBAR 

WRITE (25,603) K, IK (K) , DW (K) , RLZ , RL1 , AREA, QPZ, QP1 , WTZ , WT1 , 
& QPBAR 

140 CONTINUE 

C 

C == End inward-boundary integration. 

C 

QPAVG=QASUM/ASUM 

WRITE (25,607) ASUM, QASUM, IMON, QPAVG 
WRITE (27,705) IMON, QPAVG 
100 CONTINUE 
C 

C == End Month Loop ===========================:=========== == = ======== 

C 

CLOSE (25) 

CLOSE (27) 

END 
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APPENDIX C 

Sample Input and Output 

Examples are given for the North American study region (Figure 2.5) 

Input Files 

1. naOOl.bn -- Enumeration of the nodes defining the boundary. 

2. naOOl.sq -- The weighting fractions used in area integrals. 

Output Files 

1. naOOl.ar -- The area of the region as computed by AREA. 

2. naOOl.qef -- The components, magnitude and direction of the 

areally— averaged moisture flux vector, as computed by QEFF. 

3. naOOl.wu -- The boundary influx in the direction of (3), ( wu ) eff > as 
computed by WU. 

4. na001.txt — Segment-by-segment and net moisture flux across the 

boundary, as computed by BFLUXW. 
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naOOl .bn 


Sat Mar 30 18:10:25 1991 


1 


naOOl 


1 

54 

58 

2 

53 

58 

3 

52 

58 

4 

51 

58 

5 

50 

58 

6 

50 

57 

7 

50 

56 

8 

50 

55 

9 

50 

54 

10 

51 

54 

11 

52 

54 

12 

53 

54 

13 

54 

54 

14 

54 

55 

15 

54 

56 

16 

54 

57 

17 

54 

58 

999 

999 

999 


naOOl. sq Fri Mar 29 13:39:30 1991 1 

n a 0 C 1 

5 5 54 58 

2 4 4 4 2 

4 8 3 8 4 

4 8 8 8 4 

4 8 8 8 4 

2 4 4 4 2 
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naOOl . ar 


Fri Mar 29 15:07:27 1991 


1 


Area Computation for naOOl.sq 
0 . 1 959e+l 3 


naOOl . qef 

Fri 

Mar 29 14:08: 

: 07 1991 

Area Mean 

Q Vector for 

naO 0 1 

Month 

<Q (lam) > 

<Q (phi) > 
[ g/ s /m] 

1 <Q> 1 

1 

0 . 864e+05 

0. H2e+05 

0. 872e+05 

2 

0 . 695e+05 

0 . 417e+04 

0. 696e+05 

3 

0 . 740e+05 

0 . l84e+05 

0 . 762e+05 

4 

0 . 853e+05 

0 . 4 30e+05 

0 . 955e+05 

5 

0 . 729e+05 

0 . 360e+05 

0.81 3e+05 

6 

0 . 774e+05 

0 . 443e+05 

0 . 892e+05 

7 

0 . 765e+05 

0 . 174e+05 

0 . 784e+05 

8 

0 . 613e+05 

0 . 17 3e+05 

0 . 637e+05 

9 

0 . 807e+05 

0 . 4 59e+05 

0 . 928e+05 

10 

0 . 743e+05 

0 . 211e+05 

0 . 772e+05 

11 

0 . 854e+05 

0 . 114e+05 

0 . 862e+05 

12 

0 . 97 6e+05 

0 . 21 8e+05 

0. 100e+06 


1 


Theta 

[deg] 


7 . 4 

3.4 
14.0 
26.8 
26.3 

29.8 

12 . 8 
15.7 

29.6 
13 . 9 

7 . 6 

12 . 6 


naOOl . wu Sun May 19 11:34:45 1991 

Weighted Mean Q (Parallel) [q m'M-l) s (-1)] 


IMON <Qp> 

1 0 . 4 06e+05 

2 0 . 321e+05 

3 0 . 446e+05 

4 0 . 675e+05 

5 0 . 529e+05 

6 0 . 539e+05 

7 0 . 204e+05 

8 0 . 169e+05 

9 0 . 420e+05 

10 0 . 392e+0 5 

11 0 . 4 68e+05 

12 0 . 534e+05 


naOOl 
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na001.txt Fri Mar 29 15:12:29 1991 1 

1 naOOl 


Month 1 (jan) 


Node K 

K+l 

Flux, kg m* 
F(x) 

-1 s A -l 
F(y) 

Scaled inward 
normal vector, m 

B(x) B (y) 

Segment flux 
kg s A -l 

Qin 


1 

2 

0 . 3 42e+02 - 

0. 114e+02 

0 .278e+06 

0 . 000e+00 

0. 951e+07 

+- 

0 . 523e+0 6 

2 

3 

0 . 320e+02 - 

0. 859e+01 

0 . 278e+0 6 

0 . 000e+00 

0. 890e+07 

+- 

0 . 510e+06 

3 

4 

0 . 332e+02 - 

0 . 518e+01 

0 . 278e+06 

0 . 000e+00 

0. 922e+07 

+- 

0 . 4 98e+06 

4 

5 

0 . 433e+02 - 

0 . 160e+00 

0. 278e+06 

0. 000e+00 

0.121e+08 

+ - 

0 . 504e+06 

5 

6 

0. 635e+02 

0 . 101e+02 

0. 000e+00 

0 . 469e+06 

0 . 475e+07 

+- 

0 . 138e+06 

6 

7 

0. 958e+02 

0 . 308e+02 

0 . 000e+00 

0. 469e+06 

0 . 144e+08 

+- 

0 . 1 88e+06 

7 

8 

0 . 12 6e+03 

0 . 478e+02 

0 . 000e+00 

0 . 4 6 9e + 0 6 

0 . 224e+08 

+- 

0 . 237e+0 6 

8 

9 

0. 140e+03 

0 . 468e+02 

0 . 000e+00 

0 . 4 69e+0 6 

0 . 21 9e+08 

+ - 

0 . 240e+06 

9 

10 

0.141e+03 

0 . 382e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 392e+08 

+- 

0. 864e+06 

10 

11 

0 . 12 9e+03 

0 . 286e+02 

-0 . 278e+06 

0.000e+00 

-0 . 359e+08 

+- 

0. 904e+0 6 

11 

12 

0.111e+03 

0 . 188e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 309e+08 

+- 

0. 903e+0 6 

12 

13 

0. 921e+02 

0 . 112e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 256e+08 

+- 

0 . 8 62e+0 6 

13 

14 

0 . 748e-f 02 

0 . 3 90e+01 

0 . 000e+00 

-0 . 410e+06 

-0 . 160e+07 

+ - 

0 . 180e+0 6 

14 

15 

0 . 668e+02 - 

0 . 377e+01 

0 , 000e+00 

-0.410e+06 

0 . 155e+07 

+- 

0 . 255e+06 

15 

16 

0 . 621e+02 - 

0 . 120e+02 

0 . 000e+00 

-0 . 410e+06 

0 . 4 91e+07 

+- 

0 . 33 4e+0 6 

16 

17 

0 . 4 67e+02 - 

0 . 147e+02 

0 . 000e+00 

-0 . 410e+06 

0. 601e+07 

+- 

0 . 362e+0 6 

Area 

= 0 . 1 96e+13 m A 2 







<Qin> 

= 0.510e+01 mm/day 

+ - 0 . 542e 

-01 

Q in = 

0.116e+09 

+- 

0 . 123e+07 

<Qout> 

= 0.587e+01 mm/day 

+ - 0 . 783e 

-01 

Qout « 

0 . 133e+09 

+- 

0 . 178e+07 

<Qnet> 

= -0.770e+00 mm/day 

+- 0 . 953e 

-01 

Qnet = 

-0 . 175e+08 

+ - 

0 . 21 6e+07 


1 naOOl 


Month 2 (feb) 


Flux, kg m A -l s A -l Scaled inward Segment flux 

normal vector, m kg s A -l 


Node K 

K+l 

F(x) 

F <y) 

B(x) 

B{y) 

Qin 



1 

2 

0 . 251e+02 - 

0.10 6e+02 

0 . 278e+06 

0 . 000e+00 

0. 699e+07 

+ - 

0 . 40 9e+0 6 

2 

3 

0 . 257e+02 - 

0, 835e+01 

0 . 278e+06 

0. 000e+00 

0. 714e+07 

+ - 

0 . 417e+0 6 

3 

4 

0 . 2 74e+02 - 

■0. 474e+01 

0 . 278e+0 6 

0 . 000e+00 

0. 762e+07 

+ - 

0 . 450e+06 

4 

5 

0 . 41 0e+02 

0 . 330e+00 

0 . 278e+06 

0 . 000e+00 

0 . 114e+08 

+ - 

0 . 518e+06 

5 

6 

0. 601e+02 

0 . 788e+01 

0 . 000e+00 

0 . 4 69e+0 6 

0 . 369e+07 

+- 

0 . 1 43e+06 

6 

7 

0. 821e+02 

0 . 206e+02 

0 . 000e+00 

0. 469e+06 

0. 966e+07 

+ - 

0 . 1 60e+06 

7 

8 

0 . 110e+03 

0.310e+02 

0 . 000e+00 

0 . 4 69e+0 6 

0 . 145e+08 

+ - 

0 . 184e+06 

6 

9 

0 . 12 6e+Q3 

0 . 328e+02 

0 . 000e+00 

0 . 4 69e + 0 6 

0 . 154e+08 

+- 

0 . 188e+06 

9 

10 

0 . 120e+03 

0 . 280e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 332e+08 

+- 

0 . 8 61e+0 6 

10 

11 

0 . 104e+03 

0 . 193e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 288e+08 

+ - 

0. 836e+06 

11 

12 

0 . 880e+02 

0 . llle+02 

-0 . 278e+06 

0. 000e+00 

-0 . 245e+08 

+ - 

0 . 789e+06 

12 

13 

0 . 718e+02 

0 . 426e+01 

-0 . 278e+0 6 

0. 000e+00 

-0 . 1 99e+08 

+ - 

0 . 72 6e+0 6 

13 

14 

0 . 607e+02 - 

0. 227e+01 

0 . 000e+00 

-0 . 410e+06 

0. 930e+06 

+ - 

0 . 233e+0 6 

14 

15 

0.54 6e+02 - 

0. 969e+01 

0 . 000e+00 

-0 . 410e+06 

0 . 3 97e+07 

+ - 

0 . 31 6e + 0 6 

15 

16 

0 . 473e+02 - 

0 . 1 63e+02 

0 . 000e+00 

-0. 410e+06 

0. 668e+07 

+- 

0. 381e+06 

16 

17 

0 . 338e+02 - 

0 . 151e+02 

0 . 000e+00 

-0 . 410e+06 

0. 621e+07 

+ - 

0 . 382e*+0 6 

Area 

= 0 . 196e+13 m A 2 







<Qin> 

= 0.416e+01 mm/day 

+ - 0 . 51 7e 

-01 

Qin = 

0. 942e+08 

+ - 

0 . 117e+07 

<Qout> 

= 0.469e+01 mm/day 

+- 0 . 71 Oe 

-01 

Qout = 

0 . 106e+09 

+ - 

0 . 1 61e+07 

<Qnet> 

= -0.539e+00 mm/day 

+- 0 . 878e 

-01 

Qnet = 

-0 . 122e+08 

+ - 

0. 199e+07 
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naOOl . txt 


Fri Mar 29 15:12:29 1991 


2 


1 naOOl 

Month 3 (mar) 


Node K 


1 ' 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

Area 

<Qin> 

<Qout> 

<Qnet > 

1 naOOl 

Month 


Flux, kg m A -l s A - 1 Scaled inward 

normal vector, m 


Segment flux 
kg s A -l 


K+l 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


F(x) 


0 . 259e+02 
0 . 268e+02 
0 . 312e+02 
0 . 478e+02 
0 . 664e+02 
0 . 871e+02 
0 . 112e+03 
0 . 128e+03 
0 . 12 9e+03 
0 . 117e+03 
0.101e+03 
0 . 778e+Q2 
0 . 5 95e+02 
0 . 504e+02 
0 . 4 45e+02 
0 . 336e+02 


F(y) 


-0. 659e+01 
-0.43 6e+01 
-0 . 130e+01 
0 . 589e+01 
0 . 181e+02 
0 . 355e+02 
0 . 51 6e+02 
0 . 505e+02 
0 . 41 0e+02 
0 . 335e+02 
0 . 240e + 02 
0 . 154e+02 
0 . 102e+02 
0 . 604e+01 
-0 . 166e+01 
-0. 712e+01 


B(x) B (y) Qin 


0 . 278e+06 
0 . 278e+0 6 
0 . 27 8e+0 6 
0 . 278e+0 6 
0 . 000e+00 
0 . 000e+00 
0 . 000e+00 
0 . 000e+00 
■0 . 278e+06 
■0 . 278e+0 6 
■ 0 . 27 8e+0 6 
*0 . 278e+06 
0 . 000e+00 
0. 000e + 00 
0 . 000e+00 
0 . 000e + 00 


0 . 000e+00 
0 . OOQe+OO 
0 . 000e+00 
0 . 000e+00 
0 . 4 69e+0 6 
0 . 4 69e+06 
0 . 469e+06 
0 . 469e+06 
0 . 000e+00 
0 . 000e+00 
0 . 000e+00 
0 . 000e+00 
-0. 410e+06 
-0 . 410e+06 
-0 . 410e+06 
-0 . 410e+06 


0 . 720e+07 + - 
0 . 745e+07 + - 
0 . 867e+07 +- 
0 . 133e+08 +- 
0 . 84 9e+07 +- 
0 . 167e+08 + - 
0 . 242e+08 + - 
0 . 237e+08 +* 
-0 . 357e+08 + • 
-0 . 325e+08 +- 
-0 . 27 9e+08 + • 
-0 . 216e+08 + • 
-0 . 418e+07 + • 
-0 . 248e+07 + 
0 . 680e+06 + 
0 . 292e-f 07 + 


0. 422e+06 
0 . 443e+0 6 
0 . 4 8 0e+0 6 
0 . 547e+0 6 
0 . 1 62e+ 0 6 
0 . 193e+06 
0 . 211e+0 6 
0 . 199e+06 
0 . 798e+06 
0 . 780e+0 6 
0 . 730e+06 
0 . 661e+0 6 
0 . 175e+0 6 
0 . 20 4e+0 6 
0 . 255e + 0 6 
0 . 272e + 0 6 


= 0.19 6e+13 m A 2 

= 0.499e+01 mm/day -4— 0.481e-01 

= 0.549e+01 mm/day +- 0.667e-01 

= -0.496e+00 mm/day +- 0.822e-01 


Qin = 
Qout = 

0 . 113e+09 
0 . 12 4e+0 9 

+ - 
+ - 

0 .109e+07 
0 . 151e+07 

Qnet = 

-0 . 112e+08 

+ - 

0 . 186e+07 


4 (apr) 


Flux, kg m A -l s A -l Scaled inward 

normal vector, m 


Segment flux 
kg s A -l 


Node K 

K+l 

F (x) 

F (y) 

B (x) 

B (y ) Qin 



II 

2 

0 . 27 6e + 02 

0 . 63 8e+01 

0 . 2 7 8e+0 6 

0.0Q0e+00 0.767e+07 

+ - 

0 . 38 9e+Q 6 

2 

3 

0. 298e+02 

0 . 81 9e+01 

0 , 27 8e + 0 6 

0 . 000e+00 0 . 827e+07 

+ - 

0 . 432e+0 6 

3 

4 

0 . 33 6e+02 

0 . 106e+02 

0 . 278e+06 

0 . 000e+00 0 . 934e+07 

+ - 

0 . 4 82e + 0 6 

4 

5 

0 . 502e+02 

0 . 200e+02 

0 . 278e + 0 6 

0 . 000e+00 0 . 140e+08 

+ - 

0 . 550e + 06 

5 

6 

0 . 722e+02 

0 . 435e+02 

0 . 000e+00 

0 . 4 69e+06 0.204e+08 

+ - 

0 . 323e+ 0 6 

6 

7 

0. 921e+02 

0 . 7 4 4e + 02 

0 . 000e+00 

0 . 469e+06 0.349e+08 

+ - 

0 . 320e+0 6 

7 

8 

0 . llle+03 

0. 8 63e+02 

0 . 000e+00 

0 . 469e+06 0.405e+08 

+ - 

0 . 296e+06 

8 

9 

0 . 120e+03 

0 . 60 9e + 02 

0 . 000e+00 

0 . 469e+0 6 0.286e+08 

+ - 

0 . 252e+0 6 

9 

10 

0 . 12 9e+03 

0 . 385e+02 - 

*0 . 278e+0 6 

0 . 000e+00 -0 . 357e+08 

+ - 

0. 642e+06 

10 

11 

0 . 131e+03 

0 . 342e+02 - 

- 0 . 27 8e+0 6 

0 . 000e+00 -0 . 363e+08 

+ - 

0 . 6 68e+0 6 

11 

12 

0 . 116e+03 

0 . 25 9e+02 - 

-0 . 27 8e+ 0 6 

0 . 000e+00 -0 . 322e+08 

+ - 

0. 654e+06 

12 

13 

0. 927e+02 

0 . 179e+02 ■ 

-0 . 278e+06 

0 . 000e+00 -0 . 258e+08 

+- 

0 . 613e+0 6 

13 

14 

0 . 754e+02 

0 . 173e+02 

0 . 000e+00 

-0 . 410e+06 -0 . 709e+07 

+ - 

0 . 18 6e+0 6 

14 

15 

0 . 643e+02 

0 . 230e+02 

0 . 000e+00 

-0 . 410e+06 -0 . 943e+07 

+ - 

0 . 155e+0 6 

15 

16 

0 . 531e+02 

0 . 219e+02 

0 . 000e+00 

- 0 . 410e+06 -0 . 898e+07 

+ - 

0 . 136e+06 

16 

17 

0 .374e+02 

0 . 118e+02 

0 . 000e+00 

-0 . 410e+06 -0 . 485e+07 

+ - 

0 . 1 08e+0 6 

Area 

= 0 . 1 96e+13 m A 2 






<Qin> 

* 0 . 

722e+01 mm/day +- 0.489e 

-01 

Qin * 0.164e+09 

+ - 

0 . llle+07 

<Qout> 

= 0 . 

707e+01 mm/day +- 0.583e 

-01 

Qout * 0.160e+09 

+ - 

0 . 132e+07 

<Qnet> 

= 0. 

143e+00 mm/day +- 0.761e 

-01 

Qnet = 0 . 323e+07 

+ - 

0 . 173e+07 
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3 


1 naOOl 

Month 5 (may) 


Flux, kg m A -l s A -l Scaled inward 

normal vector, m 


Segment flux 
kg s A -l 


Node K 

K+l 

F(x) 

F(y) 

B(x) 

1 

2 

0 . 285e+02 

0 . 774e+01 

0.278e+06 

2 

3 

0 . 289e+02 

0 . 964e+01 

0 . 278e+06 

3 

4 

0 . 301e+02 

0 . 114e+02 

0 . 278e+06 

4 

5 

0. 412e+02 

0 . 220e+02 

0 . 278e+06 

5 

6 

0 . 538e+02 

0 . 458e+02 

0 . 000e+00 

6 

7 

0. 622e+02 

0. 630e+02 

0 . 000e+00 

7 

8 

0. 683e+02 

0 . 547e+02 

0. 000e+00 

8 

9 

0. 695e+02 

0 . 28 6e+02 

0 . 000e+00 

9 

10 

0 . 807e+02 

0 . 158e+02 

-0 . 278e+06 

10 

11 

0. 973e+02 

0 . 181e+02 

-0 . 278e+06 

11 

12 

0 . 103e+03 

0 . 162e+02 

-0 . 278e+06 

12 

13 

0 . 968e+02 

0 . 117e+02 

-0 . 278e+0 6 

13 

14 

0 . 890e+02 

0 . 1 4 6e+02 

0. 000e+00 

14 

15 

0. 81 0e+02 

0 . 23 4e+02 

0 . 000e+00 

15 

16 

0. 671e+02 

0 . 237e+02 

0 . 000e+00 

16 

17 

0. 431e+02 

0 . 13 6e+02 

0 . 000e+00 

Area 

= 0 . 1 96e+l 3 m A 2 



<Qin> 

■ 0.555e+01 mm/day 

+ - 0 . 3 94e 

-01 

<Qout > 

= 0.599e+01 mm/day 

+ - 0 . 41 8e 

-01 

<Qnet > 

= -Q.436e+00 mm/day 

+ - 0 . 575e 

-01 


B(y) 


Qin 


> 0.000e+00 

0 . 792e+07 

+ - 

0 . 32 6e+0 6 

> 0 . 000e+00 

0 . 802e+07 

+ - 

0 . 32 6e+0 6 

; 0 . 000e+00 

0 . 837e+07 

+ - 

0 . 337e+06 

; 0 . 000e+00 

0 . 115e+08 

+ - 

0 . 376e+06 

' 0 . 4 69e+06 

0 . 215e+08 

+ - 

0 . 283e+06 

1 0 . 4 69e+0 6 

0 . 295e+08 

+- 

0 . 292e+06 

' 0 . 4 69e+0 6 

0 . 256e+08 

+ - 

0 . 297e+06 

i 0 . 4 69e + 06 

0 . 134e+08 

+- 

0 . 282e+06 

0 . 000e+00 

-0 . 224e+08 

+- 

0 . 383e+0 6 

0 . 000e+00 

-0 . 270e+08 

+- 

0. 427e+06 

0 . 000e+00 

-0 . 285e+08 

+ - 

0. 462e+06 

0 . 000e+00 

-0 . 269e+08 

+- 

0 . 4 85e+0 6 

-0. 410e+06 

-0 . 599e+07 

+ - 

0 . 1 92e-f 0 6 

-0 . 410e+06 

-0. 957e+07 

+- 

0 . 1 8 8e+0 6 

-0. 410e+06 

-0 . 971e+07 

+ * 

0 . 177e+0 6 

-0 . 410e+06 

-0 . 5 60e+07 

+ - 

0 . 134e + 0 6 

Qin = 

0 . 12 6e+0 9 

+ - 

0 . 8 94e + 0 6 

Qout = 

0 . 136e+09 

+ - 

0 . 94 8e+0 6 

Qnet = 

-0. 98 9e+07 

+ - 

0. 130e+07 


1 naOOl 


Month 6 (jun) 


Node K 

K+l 

Flux, kg m A -l s A -l 
F(x) F(y) 

Scaled inward 
normal vector, m 

B(x) B(y) 

Segment flux 
kg s A -l 

Qin 


1 

2 

0 . 428e+02 

0 . 18 9e+02 

0 .278e+06 

0 . 000e+00 0 . 119e+08 

KSSS 

0 . 335e + 06 

2 

3 

0. 400e+02 

0 . 195e+02 

0. 278e+06 

0. 000e+00 

0 . llle+08 

+ - 

0 . 31 6e+0 6 

3 

4 

0. 367e+02 

0 . 192e+02 

0 . 278e+0 6 

0 . 000e+00 

0 . 102e+08 

+ - 

0 . 2 93e+0 6 

4 

5 

0 . 370e+02 

0 . 287e+02 

0 . 278e+06 

0. 000e+00 

0 . 103e+08 

+ - 

0 . 280e+06 

5 

6 

0 . 347e+02 

0 . 622e+02 

0 . 000e+00 

0. 469e+06 

0 . 292e+08 

+ - 

0 . 32 6e+0 6 

6 

7 

0 . 360e+02 

0. 827e+02 

0. 000e+00 

0 . 4 69e+06 

0 . 388e+08 

+ - 

0 . 360e+06 

7 

8 

0 . 450e+02 

0 , 574e+02 

0 . 000e+00 

0 . 4 6 9e + 0 6 

0 . 269e+08 

+ - 

0 . 370e + 06 

6 

9 

0. 490e+02 

0 . 173e+02 

0. 000e+00 

0. 4 69e+06 

0 . 810e+07 

+ - 

0 . 3 60e + 0 6 

9 

10 

0. 64 6e+02 

0. 407e+00 

-0 . 278e+06 

0. 000e+00 

-0 . 179e+08 

+ - 

0. 273e+06 

10 

11 

0. 942e+02 

0 . 372e+01 

-0 . 278e+06 

0 - 000e + 00 

-0 . 262e+08 

+ - 

0. 328e + 06 

11 

12 

0 . 11 9e+03 

0. 614e+01 

- 0 . 278e+0 6 

0 . 000e+00 

-0 . 329e+08 

+ - 

0 . 391e+06 

12 

13 

0 . 133e+03 

0 . 63 6e+01 

-0 . 278e+06 

0. 000e+00 

-0. 370e+08 

+ - 

0 . 4 4 9e+0 6 

13 

14 

0 . 130e+03 

0 . 936e+01 

0. 000e+00 

-0. 410e+06 

-0 . 384e+07 

+ - 

0 . 2 61e+0 6 

14 

15 

0 . 114e+03 

0 . 274e+02 

0 . 000e+00 

-0. 410e+06 

-0 . 112e+08 

+ - 

0 . 255e+0 6 

15 

16 

0. 942e+02 

0. 446e+02 

0 . 000e+00 

-0 . 410e+06 

-0 . 183e+08 

+ - 

0 . 247e+06 

1 6 

17 

0. 627e+02 

0 . 325e+02 

0 . 000e+00 

-0. 410e+06 

-0 . 133e+08 

+ - 

0 . 214e+06 

Area 

= 0 . 1 96e+13 m^2 







<Qin> 

= Q.646e+01 mm/day 

+- 0 . 413e 

-01 

Qin = 

0 . 146e+09 

+ - 

0 . 937e+0 6 

<Qout> 

= 0.709e+01 mm/day 

+- 0 . 389e 

-01 

Qout * 

0 . 161e+09 

+ - 

0 . 882e+0 6 

<Qnet> 

= -0.630e+00 mm/day 

+- 0 . 567e 

-01 

Qnet = 

-0. 143e+08 

+ - 

0 . 12 9e+07 
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1 naOOl 


Month 7 (jul) 


Node K K+l 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


Flux, leg m A 

-1 s A -l 

Scaled inw 

ard 

Segment flux 


normal vector, m 

kg s A -l 


F (x) 

F (y) 

B (x) 

B (y) 

Qin 


0. 494e+02 

0 . 642e+01 

0 . 278e+06 

0 . 000e+00 

0 . 137e+08 

+ - 

0 . 385e+02 

0 . 777e+01 

0 . 278e+06 

0 . 000e+00 

0 . 107e+08 

+- 

0 . 220e+02 

0 . 782e+01 

0 . 278e+0 6 

0. 000e+00 

0 . 613e+07 

+- 

-0 . 750e+01 

0 . 141e+02 

0 . 278e+0 6 

0 . 000e+00 

-0.208e+07 

+- 

-0 . 204e+02 

0. 420e+02 

0 . 000e+00 

0. 469e+06 

0 . 197e+08 

+ ~ 

0 . 475e+01 

0 . 548e+02 

0 . 000e+00 

0 . 469«+06 

0 . 257e+08 

+ - 

0 . 427e+02 

0 . 2 91e+02 

0 . 000e+00 

0 . 469e+06 

0 . 136e+08 

+- 

0 . 703e+02 

0 . 71 4e+01 

0 . 000e+00 

0. 469e+06 

0 . 335e+07 

+ - 


0 . 903e-f 02 
0.111e+03 
0 . 127e+03 
0 . 138e+03 
0 . 138e+03 
0 . 12 9e+03 
0 . 109e+03 
0 . 737e+02 


-0 . 28 9e+01 
-0. 92 8e+01 
-0 . 118e+02 
-0.132e+02 
-0. 124e+02 
-0 . 303a+01 
0 . 101e+02 
0 . 101e+02 


-0. 278e+06 
-0 . 278e+06 
-0 . 278e+06 
-0.278e+06 
0. 000e+00 
0 . 000e+00 
0 . 000e+00 
0 . 000e+00 


0 . 000e+00 
0 . QQOe+OO 
0 . 000e+00 
0. 000e+00 
-0 . 410e+06 
**0 . 4l0e+06 
-0 . 410e+06 
-0 . 410e+06 


-0 . 309e+08 
-0. 352e+08 
-0 . 382e+08 
0 . 510e+07 
0 . 124e+07 
-0 . 413e+07 
-0. 412e+07 


0 .330e+06 
0.24 9e+0 6 
0 . 169e+06 
0 . 175e+06 
0 . 375e+06 
0. 407e+06 
0 . 427e+0 6 
0 . 4 17e+0 6 
0 . 27 9e+0 6 
0 . 325e+0 6 
0 . 378e+06 
0 . 42 6e+0 6 
0 . 32 9e+0 6 
0 . 351e+0 6 
0 . 32 6e+06 
0 . 227e+0 6 


Area = 

<Qin> = 
<Qout> = 

<Qnet> = 

1 naOOl 


0 . 1 96e+13 m A 2 

0 . 4 3 8e+01 mm/day +- 0.461e-01 
0 . 61 6e+01 mm/day + - 0.368e-01 

-0.178e+01 mm/day + - 0.590e-01 


Qin = 0 . 993e+08 + - 0.105e+07 
Qout = 0 . 140e+0 9 + - 0.834e+06 


Qnet 


= -0.404e+08 +- 0.134e+07 


Month 8 (aug) 


Node K K+l 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


Flux, kg m 

A -1 s A -l 

Scaled inward 

Segment flux 


normal vector, m 

kg s A -l 


F (x) 

F(y) 

B (x) 

B (y) 

Qin 

EXCS 

0. 368e+02 

0 . 32 6e+01 

0 . 27 8e+0 6 

0 . 000e+00 

0 . 102e+08 

+ ~ 

0 . 290e+02 

0 . 491e+01 

0 . 278e+06 

0 . 000e+00 

0 . 808e+07 

+ - 

0 . 1 83e+02 

0 . 541e+01 

0 . 27 8e+0 6 

0 . 000e+00 

0 . 507e+07 

+ “ 

-0 . 235e+01 

0 . 929e+01 

0 . 278e+06 

0 . 000e+00 

-0 . 653e+06 

+ ~ 

-0 . 915e+01 

0 . 396e+02 

0 . OOOe+OQ 

0 . 469e+06 

0 . 18 6e+08 

+ ~ 

0. 645e+01 

0 . 566e+02 

0 . 000e+00 

0. 469e+06 

0 . 265e+08 

+ - 

0 . 274e+02 

0 . 287e+02 

0. 000e+00 

0 . 4 69e+06 

0 . 135e+08 

+ ~ 

0 . 50 9e+02 

0 . 5 93e+01 

0 . 000e+00 

0. 469e+06 

0 . 278e+07 

+ “ 


0.751e+02 
0. 990e+02 
0 . 118e+03 
0 . 131e+03 
0 .124e+03 
0 . 105e+03 
0 . 8 68e+02 
0 . 572e+02 


*0 . 150e+01 
-0. 724e+01 
'0 . 108e+02 
-0 . 13 9e+02 
-0 . 1 60e+02 
-0.575e+01 
0 . 100e+02 
0 . 881e+01 


-0 . 278e+06 
-0 . 278e+06 
-0 . 278e+06 
-0 . 278e+06 
0. 000e+00 
0 . 000e+00 
0 . 000e+00 
0. 000e+00 


0 . 000e+00 
0 . 000e+00 
0 . 000e+00 
0 . 000e+00 
-0 . 410«+0 6 
-0 . 410e+06 
-0.410e+06 
-0 . 410e+06 


0 . 209e+08 
-0.275e+08 
-0 . 327e+08 
-0 . 364e+08 
0 . 658e+07 
0 . 23 6e+07 
-0. 410e+07 
-0 . 3 61e+07 


Area = 0.196e+13 m A 2 

<Qin> = 0 . 413e+01 mm/< 

<Qout> = 0 . 555e+01 mm/< 

<Qnet> » -0 . 142e+01 mm/< 


25 4e+0 6 
1 99e+0 6 
151e + 0 6 
175e+06 
3 71e+0 6 
3 8 9e+0 6 
, 407e+0 6 
. 411e+0 6 
. 27 9e+06 
. 323e+06 
. 37 9e+06 
. 43 4e+0 6 
. 310e+06 
. 314e+06 
. 292e+0 6 
. 21 4e+0 6 


+ - 
+ - 

0. 429e-01 
0 . 3 62e- 01 

Qin * 
Qout * 

0 . 937e+08 
0 .126e+09 

+ - 
+ - 

0 

0 

+ - 

0.562e-01 

Qnet - 

-0 . 322e+08 

+ - 

0 


- 0 . 127e+07 
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1 naOOl 


Month 9 (sep) 


Node K 

K+l 

Flux, kg m / 
F(x) 

V -1 s A -l 
F(y> 

Scaled inward 
normal vector, m 

B(x) B (y) 

Segment flux 
kg s A -l 

Qin 


1 

2 

0. 487e+02 

0 . 1.05 •+ 02 

0 . 278e+06 

0 . 000e+00 

0 . 135e+08 

II i 

u + 

u 

0 . 359e+0 6 

2 

3 

0. 443e+02 

0 . 133e+02 

0 . 278e+0 6 

0. 000e+00 

0 . 123e+08 

+- 

0 . 322e+0 6 

3 

4 

0 . 372e+02 

0. 146e+02 

0. 278e+06 

0 . 000e+00 

0 . 103e+08 

+- 

0 . 278e+0 6 

4 

5 

0 . 295e+02 

0 . 223e+02 

0 . 278e+06 

0. 000e+00 

0. 819e+07 

+- 

0 . 251e+0 6 

5 

6 

0 . 2 62e+02 

0 . 438e+02 

0 . 000e+00 

0. 469e+06 

0 . 205e+08 

+- 

0 . 332e+0 6 

6 

7 

0 . 250e+02 

0 . 578e+02 

0 . 000e+00 

0 . 4 69e+06 

0 . 271e+08 

+- 

0 . 35 6e+0 6 

7 

8 

0 . 199e+02 

0 . 489e+02 

0. 000e+00 

0 . 4 69e+0 6 

0 . 229e*f 08 


0 . 375e+06 

e 

9 

0 . 151e+Q2 

0 . 272e+02 

0 . 000e+00 

0 . 4 69e+06 

0 . 128e+08 

+- 

0 . 367e+06 

9 

10 

0 . 327e+02 

0 . 173e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 909e+07 

+- 

0 . 275e+06 

10 

11 

0. 746e+Q2 

0 . 231e+02 

-0 . 278e+06 

0. 000e+00 

-0 . 207e+08 

+- 

0 . 334e+06 

11 

12 

0.111e+03 

0 . 271e+02 

-0 . 278e+06 

0. 000e+00 

-0 . 308e+08 

+- 

0. 403e+0 6 

12 

13 

0 . 133e+03 

0 . 283e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 368e+08 

+- 

0. 467e+06 

13 

14 

0 . 145e+03 

0 . 350e+02 

0 . 000e+00 

-0 . 410e+0 6 

-0 . 143e+08 

+- 

0 . 222e-f 0 6 

14 

15 

0 . 1 40e+03 

0. 459e+02 

0 . 000e+00 

-0. 410e+06 

-0 . 188e+08 

+- 

0 . 22 7e+0 6 

15 

16 

0 . 115e+03 

0. 433e+02 

0. 000e+00 

-0. 410e+06 

-0 . 177e+08 

+- 

0 . 21 9e+0 6 

1 € 

17 

0 . 749e+02 

0 . 225e+02 

0 . 000e+00 

-0 . 410e+06 

-0. 924e+07 

+ - 

0 . 1 63e+06 

Area 

= 0 . 1 96e+13 m A 2 







<Qin> 

= 0.563e+01 mm/day 

+ - 0 . 415e 

-01 

Qin = 

0 . 128e+09 

+- 

0. 941e+06 

<Qout> 

= 0.695e+01 mm/day 

+ - 0 . 380e 

-01 

Qout = 

0 . 158e+09 

+- 

0 . 8 62e+0 6 

<Qnet> 

= -0.132e+01 mm/day 

+- 0 . 563e 

-01 

Qnet = 

-0. 299e+08 

+ - 

0 . 12 8e+07 

1 naOOl 










Month 10 (oct) 



Flux, kg m~-l s A -l 

Scaled inward 

Segment flux 





normal vector, m 

kg s A -l 



Node K 

K+l F (x) 

F(y) 

B (x) 

B (y) 

Qin 



1 

2 0 . 373e+02 - 

•0 . 383e+01 

0 .278e+06 

0. 000e+00 

0 . 104e+08 

+ - 

0 . 352e+06 

2 

3 0 . 350e+02 

0 . 350e+00 

0 . 278e+0 6 

0. 000e+00 

0 . 973e+07 

+ - 

0 . 32 4e+0 6 

3 

4 0 . 325e+02 

0. 487e+01 

0. 278e+06 

0. 000e+00 

0. 902e+07 

+ - 

0 . 3 02e+0 6 

4 

5 0 . 354e+02 

0 . 154e+02 

0 . 278e+06 

0 . 000e+00 

0. 984e+07 

+ - 

0 . 303e+06 

5 

6 0 . 463e+02 

0 . 271e+02 

0 . 000e+00 

0. 469e+06 

0 . 127e+08 

+- 

0 . 224e+06 

6 

7 0 . 541e+02 

0 . 269e+02 

0 . 000e+00 

0. 469e+06 

0 . 12 6e+08 

+ - 

0 . 24 7e+0 6 

7 

8 0 . 504e+02 

0. 174e+02 

0 . 000e+00 

0. 4 69e + 0 6 

0. 816e+07 

+ - 

0 . 262e+0 6 

8 

9 0 . 435e+02 

0. 109e+02 

0 . 000e+00 

0. 469e+06 

0 . 50 9e+07 

+ - 

0 . 257e+0 6 

9 

10 0 . 538e+02 

0 . 1 41e+02 

-0 . 278e+0 6 

0 . 000e+00 

-0 . 149e+08 

+ - 

0 . 380e+0 6 

10 

11 0 . 763e+02 

0 . 220e+02 

-0 .278e+06 

0 . 000e+00 

-0 . 212e+08 

+ - 

0 . 438e+0 6 

11 

12 0 . 917e+02 

0 . 285e+02 

-0 . 278e+06 

0 . 000e+00 

-0 . 255e+08 

+- 

0 . 4 93e+06 

12 

13 0 . 102e+03 

0 . 326e+02 

-0 . 278e+06 

0. 000e+00 

-0 . 283e+08 

+ - 

0 . 536e+06 

13 

14 0 . 106e+Q3 

0 . 314e+02 

0. 000e+00 

-0 . 410e + 0 6 

-0 . 129e+08 

+ - 

0 . 168e+06 

14 

15 0 . 101e+03 

0 . 247e+02 

0 . 000e+00 

-0 . 410e+06 

-0 . 101e+08 

+- 

0 . 1 86e+06 

15 

16 0 . 873e+02 

0 . 12 9e+02 

0 . 000e+00 

-0. 410e+06 

-0 . 531e+07 

+ - 

0 . 212e+0 6 

1 6 

17 0 . 577e+02 - 

0 . 195e+00 

0 . 000e+00 

-0. 410e+06 

0 . 799e+05 

+ - 

0 . 202e+0 6 

Area 

* 0 . 1 96e+13 m A 2 







<Qin> 

- 0.342e+01 mm/day 

+- 0 . 369e 

-01 

Qin = 

0. 776e+08 

+ - 

0 . 83 6e+06 

<Qout> 

= 0,522e+01 mm/day 

+- 0 . 435e 

-01 

Qout = 

0 . 118e+09 

+ - 

0. 987e+0 6 

<Qnet> 

* -0.179e+01 mm/day 

+- 0 . 571e 

-01 

Qnet « 

-0. 407e+08 

+ - 

0 . 12 9e+0 7 
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6 


1 naOOl 


Month 11 (nov) 


Node K 

K+l 

Flux, kg m 
F (x) 

A -1 s A -l 
F(y) 

Scaled inward 
normal vector, m 

B (x ) B(y) 

Segment flux 
kg s A -l 

Qin 


1 

2 

0 . 383e+02 

-0. 840e+01 

0 . 27 8e+0 6 

0. OOOo+OO 

0.107e+08 

+ - 

0 . 4 45e+0 6 

2 

3 

0 . 373e+Q2 

-0 . 596e+01 

0. 2780+06 

O.OOOo+OO 

0.104o+08 

+ - 

0 . 443e+0 6 

3 

4 

0 . 391e+02 

-0. 2390+01 

0 . 278e+06 

0 . OOOo+OO 

0.109o+08 

+ - 

0. 451e+06 

4 

5 

0 . 532e+02 

0 . 510e+01 

0. 2780+06 

O.OOOo+OO 

0. 1480+08 

+ - 

0. 484e+06 

5 

6 

0. 731e+02 

0. 156e+02 

O.OOOo+OO 

0. 469o+06 

0. 7340+07 

+ - 

0 . 159e+06 

6 

7 

0. 9370+02 

0. 275o+02 

‘O.OOOo+OO 

0. 4690+06 

0. 1290+08 

+ - 

0. 1830+06 

7 

8 

0.109O+03 

0. 3420+02 

0 . OOOo+OO 

0 . 4690+06 

0.160O+08 

+ - 

0 . 200O+06 

8 

9 

0 . 110e+03 

0 . 2860+02 

0 . 000e+00 

0 . 469o+06 

0. 1340+08 

+ - 

0 . 193e+06 

9 

10 

0 . llle+03 

0 . 22 8e+02 

-0 . 278e+0 6 

O.OOOo+OO 

-0.309e+08 

+ - 

0 . 6590+06 

10 

11 

0 . 113e+03 

0. 213e+02 

-0 . 278e+06 

O.OOOo+OO 

-0. 313e+08 

+ - 

0 . 68 0e+0 6 

11 

12 

0.107e+03 

0 . 1 81e+02 

-0 . 278e+0 6 

0 . OOOo+OO 

-0 .297e+08 

+- 

0. 672e + 0 6 

12 

13 

0 . 978e+02 

0 . 1 43e+02 

-0 . 27 8e+0 6 

0 . 000e+00 

-0 . 272e+08 

+ - 

0 . 65 0e+0 6 

13 

14 

0 . 864e+02 

0. 736e+01 

0 .000e+00 

-0 . 410e+06 

-0.301e+07 

+ - 

0. 166e+06 

14 

15 

0 . 77 9e+02 

-0 . 320e+00 

0 . OOOe + OQ 

-0 . 410O+06 

0 . 1310+06 

+ - 

0. 207e + 06 

15 

16 

0 . 693e+02 

-0. 673e+01 

0 . 000e + 00 

-0 . 410e + 06 

0 . 27 6e+07 

+ - 

0 . 250e + 0 6 

16 

17 

0 . 515e+ 02 

-0. 982e+01 

0 . 000e+00 

-0 . 410e+06 

0 . 403e+07 

+ - 

0 . 2 57e+0 6 


Area = 0.196e+13 2 

<Qin> = 0.45 6e+01 mm/day +- 0.471e-01 Qin = 0.103e+09 + - 0.107e+07 

<Qout> = 0 . 538e+01 mm/day +- 0.591e-01 Qout = 0.122e+09 + - 0.134e+07 

<Qnet> - -0 . 83 Qe + 00 mm/day + - 0.756e-01 Qnet = -0.188o+08 +- 0.171e+07 

1 naOOl 


Month 12 (dec) 


Node K 

K+l 

Flux, kg m 
F(x) 

A -1 s A -l 

F (y) 

Scaled inward 
normal vector, m 

B (x) B(y) 

Segment flux 
kg s A -l 

Qin 


1 

2 

0.380e+02 

-0. 4490+01 

0 . 278e+0 6 

0 . OOOo+OO 

0 . IO60+O8 

+ - 

0 . 530e+0 6 

2 

3 

0 . 378e+02 

-0. 149e+01 

0 . 278e + 0 6 

0 . OOOo+OO 

0 . 105o+08 

+- 

0 . 542o+06 

3 

4 

0. 3860+02 

0. 2150+01 

0 . 278e+06 

0 . OOOo+OO 

0.107e+08 

+ - 

0 . 5 48e+06 

4 

5 

0.520O+02 

0. 966e+01 

0 . 278e+06 

0 . OOOo+OO 

0 . 1440+08 

+ - 

0.56 9e+0 6 

5 

6 

0 . 7 61e+02 

0 . 21 9e+02 

0 . OOOo+OO 

0 . 4 69e+06 

0.103O+08 

+- 

0 . 193e+06 

6 

7 

0 . 108e+03 

0. 450O+02 

0 . 000e+00 

0 . 469e+06 

0.211O+08 

+- 

0 . 248e+06 

7 

8 

0 . 136e+03 

0. 633e+02 

0 . 000e+00 

0. 4690+06 

0 . 297e+08 

+ - 

0 . 27 9e+06 

6 

9 

0 . 1 42e+03 

0 . 557e+02 

0 . 000e+00 

0 . 469e+06 

0 . 261e+08 

+- 

0 . 2 4 4e+0 6 

9 

10 

0 . 1450+03 

0. 425o+02 

-0 . 278e + 0 6 

O.OOOo+OO 

-0. 403e+08 

+ - 

0 . 807e+0 6 

10 

11 

0. 144e+03 

0 . 342e+02 

-0 . 278e+0 6 

0 . OOOo+OO 

-0. 399e+08 

+ - 

0 . 87 6e+0 6 

11 

12 

0 . 12 6e+03 

0 . 24 6e + 02 

-0 . 278e+06 

0 . OOOo+OO 

-0.350O+08 

+ - 

0 . 8 8 6e + 0 6 

12 

13 

0 . 104e+03 

0 . 1 60e+Q2 

-0 . 278e + 06 

0 . 000e+00 

-0 . 288e+08 

+ - 

0 . 83 6e + 0 6 

13 

14 

0. 873e+02 

0. 874e + 01 

0 . 000e+00 

-0 . 410e+06 

-0. 358o+07 

+ - 

0 . 127e+06 

14 

15 

0 . 766e+02 

0 . 2 61e+01 

0 . 000e+00 

-0 . 410o+06 

-0. 107O+07 

+- 

0 . 1 4 4e+0 6 

15 

16 

0. 65 4e + 02 

-0 . 354e+01 

0 . 000e+00 

-0. 410O+06 

0. 1450+07 

+ - 

0 . 1 72e+Q 6 

16 

17 

0 . 48 8e+02 

-0 . 647e+01 

0 . 000e+00 

-0. 4 1 0e+06 

0 . 265e+07 

+ - 

0 . 1 77e+0 6 


Area = 0.196e+13 m A 2 

<Qin> ■ 0 . 606e+01 mm/day +- 0.539e-01 Qin = 0.138e+09 + - 0.122e+07 

<Qout> = 0 . 656e+01 mm/day + - 0.756e-01 Qout = 0.149e+09 + - 0.171e+07 

<Qnet> * -0.4 92e+00 mm/day + - 0.929e-01 Qnet = -0.111e-f08 +- 0.211e+07 
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